wo 01/20477 



'7 



PCT/DEOO/02930 



Computerbasiertes Verfahren zur automatischen Aufbereitung 
von Daten biomagnetischer Felder, 
insbesondere von magnetokardiographischen Oaten 

5 Die Erfindung betrifft ein computerbasiertes Verfahren zur autonnatischen Aufbereitung 
von Daten biomagnetischer Felder, insbesondere von magnetokardiographischen Da- 
ten. Insbesondere betrifft die Erfindung ein Verfahren, das es eriaubt, da& mittels eines 
Oder mehrerer Detektoren, insbesondere mittels eines SQUID-Detektors in einer Ebe- 
ne bzw. an mehreren Punkten einer Ebene aufgenommene biomagnetische Felder 

10 automatisch z.B. in Form von Magnetfeldkarten zu visualisieren und die in dem Feld 
enthaltenen Informationen automatisch zu evaluieren und zu klassifizieren. 

Da sich biomagnetischen Felder - anders als die elektrischen Felder - 
nicht-leitungsgebunden ausbreiten und weitgehend unverfalscht auderhalb eines bio- 
logischen K6rpers, also zum Beispiel in einer Ebene uber dem menschlichen Herzen 

15 nicht-invasiv gemessen werden konnen, stellt die Erfassung biomagnetischer Felder 
mittels neuer, hoch-sensitiver ErfassungsgerSte, wie zum Beispiel dem in der interna- 
tionalen Patentanmeldung PCT/DEOO/02472 beschriebenen Magnetographen, der 
vorteilhaft in einer nicht-abgeschirmten Umgebung eingesetzt werden kann, eine 
Technik dar, die in den kommenden Jahren breite medizinische Anwendung finden 

20 wird. Insbesondere werden Vorsorgeuntersuchungen mdglich, bei denen in einfacher 
und schneller Weise weite BevdIkerungskreise zum Beispiel auf ihr Herzinfarktrisiko 
untersucht werden kfinnen. 

Die Entwicklung hoch-sensitiver Magnetographen, insbesondere solcher, die in 
nicht-abgeschirmten Umgebungen eingesetzt werden kdnnen. erfordert neue Verfah- 

25 ren zur Aufbereitung der erfaflten Daten, da die bislang bekannten Verfahren 3u&erst 
zeitaufWendig sind und nur von wenigen Spezialisten ausgefOhrt werden k6nnen. Die 
Daten sollen weitestgehend automatisch so aufbereitet werden, da& die in den erfaU- 
ten Magnetfeldern enthaltenen diagnotisch relevanten Informationen von einem Arzt 
leicht erkannt und bewertet werden kOnnen. Bei einer bevorzugten Ausfuhmngsform 

30 soli das Verfahren die in den Daten enthaltenen diagnostisch relevanten Informationen 
mit gespeicherten Infomiationen vergleichen und automatisch klassifizieren, um so den 
Arzt bei seiner Diagnose zu unterstutzen. Bestimmte Infomnationen sollen z.B. in Form 
von Stromlinienkarten visualisiert werden. 

Sollen biomagnetische Stroma, insbesondere die Strome im Herzen, die in der 

.35 klinischen Diagnostik enwiesenermalien eine besonders hohe Relevanz Ijesitzen, dar- 
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gestellt werden. muli ein sogenanntes "inverses Problem" geldst werden, das heilit, 
da& aus den erfadten Magnetfeldern, die von den biomagnetischen StrSmen hervorge- 
rufen werden, die die Felder hervorrufenden StrSme berechnet werden mOssen. 
In Anbetracht der KomplexitSt eines inversen Problems bei der Lokalisation 
5 magnetischer Feldquellen werden verschiedene vereinfachende Annahmen bei der 

L6sung gemacht. Eine dieser Annahmen ist beispielsweise, da& alle Quellen in dersel- 
ben Ebene positioniert und in einer bestimmten Distanz von einer Melioberfiache an- 
geordnet sind. Die erzielten Resultate davon ab. welches physikallsche Modell und 
entsprechend welches mathematische Modell man als Modell der magnetischen Feld- 
10 quellen wahlt. Falls man einen fiachen Quellenbereich wahit, wird ublichenveise die 
Verteilung der Stromdichten auf einer OberflSche als Modell der Quellen gewShlt, die 
die erfaflte Vektorgrolie darstellt. Nach Losung eines inversen Problems erhSIt man 
eine Karte von gerichteten Pfeilen, deren Grolie proportional zum Betrag einer Strom- 
dichte ist. 

15 Das Abbilden von Str6men in einer solchen Weise bedingt bestimmte Schwie- 

rigkeiten. Auf Seiten des berechneten Bereiches, wo der Modulo-Wert der Stromdich- 
ten mehrere Male kieiner ist als das Maximum, ist es nicht ohne weiteres mOglich, auf 
der Karte von Pfeilen die Gr5(ie und Richtung der Pfeile abzuschatzen. Die Quellen mit 
kieiner Intensitat fallen so aus der Untersuchung. 

20 Ein Verfahren zur automatischen Visualisieren von Herzstromen durch Stromli- 

nien in der Ebene wird in dem Aufsatz von R, Killman. G.G. Jaros, P. Wach, R. Grau- 
mann. W. Moshage. M. Renhard, P.H. Fleischmann: Localisation of Myocardial 
Ischaemia from the Magnetocardiogramm using Current Density Reconstruction Me- 
thod: Computer Simulation Study. Biological Engineering & Computing, September, 

25 1 995, Seiten 643 - 651 , beschrieben. GemSIl dem in dem genannten Aufsatz be- 

schriebenen Verfahren erhait man als Resultat der Ldsung des inversen Problems eine 
der genannten Karten von orientierten Pfeilen. deren GrSlie proportional zum Modulus 
der Stromdichte ist. 

Das aus dem genannten Aufsatz bekannte Verfahren zum Abbilden von Herz- 

30 stromen durch Stromlinien in der Ebene hat jedoch den Nachteil, dad in Bereichen, in 
denen die Stromdichten modulu 3 und mehr geringer sind als das Maximum, eine vi- 
suelle Beurteilung des mit dem Verfahren erzeugten Musters von Pfeilen hinsichtlich 
ihrer Gr6(le und Richtung nicht mehr ohne weiteres mOglich ist. In der Praxis hat sich 
gezeigt, dali im Hinblick auf eine vernunftige Fehlerabschatzung lediglich ein bis zwel 

35 besonders intensive Bereiche eine magnetischen Feldes visualisiert werden k5nnen. 
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Die Obrigen. weniger intensiven Bereiche konnen nicht sichtbar wahrgenommen wer- 
den. 

Davon ausgehend liegt der Erfindung zum einen die Aufgabe zugrunde. ein 
Verfaliren zum Aufbereiten der Daten des erfaliten biomagnetischen Feldes anzuge- 
ben, mittels welchem die das biomagnetische Feld hervorrufenden Strfime durch 
Stromlinlen in der Ebene visualisiert werden kdnnen und so einer visuellen Auswertung 
durch geschultes Personal Oder einen Arzt in besonders einfacher Weise zuganglich 
werden. 

Die Aufgabe wird gel6st yon einem Verfahren mit den IVlerkmalen des An- 
spruchs 1 . 

Als Modell der biologischen Strome wird erfindungsgemafi die Oberflachendich- 
te einer doppelten Schicht von magnetischen Ladungen (einfache Schicht magneti- 
scher Dipole) oder, was aquivalent hierzu ist. eine Funktlon der Strome gewahlt. Eine 
bemerkenswerte Eigenschaftdieser Funktion ist, da(i die Projektion Ihrer Niveaulinien 
auf eine berechnete Ebene Stromlinien darstellt. Deshalb erhait man durch Losen einer 
Integralgleichung betreffend die Dichte einer doppelten Schicht von magnetischen La- 
dungen eine Karte der Stromlinien, wobei die rechten Glieder der Gleichung die Vertei- 
lung des z.B. in einer Ebene nahe des Thorax gemessenen magnetischen Feldes 
(Komponenten des magnetischen Induktionsvektors oder dessen Ableitungen) wieder- 
geben. 

Erflndungsgemafi wird angenommen, dad die tatsdchlich voluminSsen Quellen 
des Feldes in einer (berechneten) Ebene angeordnet sind, die parallel zu der gemes* 
senen Ebene verlSuft. GemaH einem akzeptierten physikalischen und mathematischen 
Modell der Quellen eines magnetischen Feldes wird eine Integralgleichung aufgestellt 
und geiast, deren rechtes Glied das gemessene magnetische Feld darstellt und deren 
Ldsung die Verteilung seiner Quellen im untersuchten KOrper, also z.B. Im Herzen Ist. 
Falls diese Quellen den Charakter separater magnetischer Dipole besitzen, deren Ach- 
se senkrecht zu einer Ebene gefQhrt ist, werden sie als magnetische Blatter interpret 
tiert und in einer unter Verwendung des Verfahrens auf einer Ausgabeeinheit, z.B. ei- 
nem Bildschirm oder einem Dmcker ausgegebenen, vorzugsweise farbigen Karte. se- 
hen sie aus wie separate Stromwirbet. 

Die wesentliche Neuerung liegt in der Behandlung des physikalischen und ma- 
thematischen Modells der Quellen eines elementaren magnetischen Feldes in einer 
Ebene als Funktion des Stroms und seine VenA/endung in einer Integralgleichung als 
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gesuchte Unbekannte. Die nach Ldsung dieser Gleichung erhaltenen Niveaulinien der 
Funktion (die Linien konstanter Werte) reprSsentieren die Stromlinien.. 

Die Karten der Stromlinien geben im Vergleich zu den Karten der Stromdiclite- 
vektoren (Pfeile) eine mehr visuelle Representation der Ausbreitung eines Stroms in 
5 einer berechneten Ebene. Auf einer Karte von Stromlinien sind die Quellen sowohl 
grolierals auch kleiner Intensitat prazise abgebildet. Die Stromwirbel, die sich nahe 
berechneten Bereichsorten befinden, werden unabhangig von ihrer Gr5(ie und Intensi- 
tat abgebildet. Deshalb macht es das Verfahren moglich, nicht nur ubiiche Muster der 
Biostrdme zum Beispiel im Herzen zu sehen, sondern auch Details dieses Stromflus- 

10 ses zu unterscheiden. 

Die Resultate einer Losung des inversen Problems durch das genannte Verfah- 
ren eriaubt es, Karten von Stromdichtevektoren aufzubauen. wahrend die Resultate, 
die mit bislang bekannten Verfahren erhalten werden, nicht dazu verwendet werden 
konnen, Stromlinienkarten aufzubauen. 

15 Die Stromfunktion ist ein Skalar, und die Stromdichte ist ein Vektor. der in der 

Ebene zwei Projektionen hat. Deshalb ist fUr dasselbe Ebenennetz die Menge ben6tig- 
ter Werte einer Stromfunktion zweimal kleiner als die Menge der benOtigten Werte bei 
der Betrachtung der Projektionen von Stromdichtevektoren. Dementsprechend ist auch 
die Ordnung des zu Idsenden linearen algebraischen Gleichungssystems zweimal 

20 niedriger. 

Nach L6sung des algebraischen Systems und Differenziation einer Stromfunkti- 
on ist es m5glich, die Verteilung der Stromdichtevektoren wie in der traditionellen For- 
mulierung zu finden. 

Die Abbildung der Strome in der Ebene durch Stromlinien eriaubt es auf, die 

25 Genauigkeit der Ldsung wesentlich zu verbessem, so dafi es m5glich wird, auch eher 
kleine Details einschlie&lich der Stromwirbel kleinen Durchmessers zu visualisieren. 

Handelt es sich bei dem erfaliten Daten um magnetokardiogrphische 
(MCG)-Daten, so sollen elektrische Ereignisse im Herzen evaliert und Ruckschlusse 
auf den physiologischen Status des Herzens gezogen werden. Die fur diesen Zweck 

30 benotigten Informationen sind in dem oberhalb der Thorax-Oberflache aufgezeichneten 
magnetischen Feld des Herzens enthalten, jeodch auf Gnjnd ihrer "kryptischen" Form 
nicht fur eine direkte Interpretation geeignet. Eine besondere Aufgabe der Erfindung ist 
daher die Aufbereitung solcher sich wShrend des Kardiozyklus andernden magneto- 
kardiographischer Daten derart, dali ein l\4ediziner aus ihnen direkt Aufschluli uber den 

35 physiologischen Zustand des untersuchten Herzens erhalt. Es sollte dazu ein haltbares 
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und klares System von diagnostischen Kriterien entwickelt werden, vergleichbar dem 
BUS der Elektrokardiologie (EK6) bekannten Standardkriteriensatz. 

Dazu wird ein computerbasiertes Verfahren zur automatischen Aufbereitung 
von Daten biomagnetischer Felder, insbesondere von magnetokardiographischen Da- 
5 ten, vorgeschlagen, bei welchem zumindest ausgewShlte Daten nach vorgegebenen 
Kriterien automatisch mit wenigstens einem vordefinierten Nomnalwert vergtichen wer- 
den und dad bei Abweichung der Daten von dem Normalwert um ein vorgegebenes 
IVlali ein auf einer Ausgabeeinrichtung ausgebbares die Abweichung signalisierendes 
Signal erzeugt wird. Das Signal kann akustisch, insbesondere aber optisch ausgege- 

10 ben werden, z.B. durch Darstellung eines aus den Daten ermittelten Wertes in rot, 
wahrend Normalwerte in schwarz oder grun dargestellt werden. 

ErfindungsgemaS durchlaufen die erfaliten Daten dazu automatisch unter- 
schiedliche Analysestufen, wobei die Komplexitat der Analyse auf jeder Stufe zunimmt. 
Die einzelnen Stufen sind in Tabelle 1 zusammengaUt. Im einzelnen werden folgende 

15 Analysestufen durchlaufen: 

1 . Amplituden-Zeit-Analyse verschiedener MCG-Kurven. 

2. Analyse hoch-aufldsender MCG. 

3. Analyse summierter MCG und Vector MCG (VMCG). 

4. Qualitative Analyse der magnetischen Feld-Verteilungskarten. 
20 5. Quantitative Analyse der magnetischen Feld-Verteilungskarten. 

6. Analyse der effektiven Dipol-Charakteristlken. 

7. Analyse der zweidimensionalen Ladungsverteilung. 

8. Analyse der dreidimensionalen Ladungsverteilung. 

Die erste Analysestufe ist 3hnlich dem Verfahren bei der morphologischen 
25 Analyse Qbllcher EKG, zumal MCG-Kurven wie EKG-Kurven aussehen und beide die- 
selbe Nomenklatur von Weilen und Intervallen - PQRST - besitzen. Auf dieser Analy- 
sestufe kann berelts myocardial IschSmie erkannt werden. 

Bel der zweiten Analysestufe, der spektrotemporale Analyse (d.h. der Bestim- 
mung der relativen Energie eines Kardiosignal-Spektrums fur verschiedene Frequenz- 
30 bander und Bestimmen der spektralen Variabilitat) und Zeitbereich-Analyse (vor allem 
QRS-Dauer) des Signals) wird in verschiedenen Mefipunkten ein Durchschnitt des 
MCG gebildet. Ublicherweise ist der Zweck einer solchen Analyse die Bestimmung der 
Homogenitat der ventrikularen Depolarisation und Verwendung dieser Daten zur Eva- 
luation des Risikos des Auftretens von Arhythmien. Ein solcher Ansatz wird auch an- 
35 gewandt um das Risiko der AbstoBung von Transplantaten abzuschatzen. 
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Auf der dritten Analysestufe wird uber alie MeBpunkten summiert. Ein solcher 
Ansatz gibt eine mehr generalisierte RepSsentation uber einige Eigenschaften der 
myocardialen Anregung. Vorteilhaft ist es insbesondere. die Felder unter der P-Welle 
und dem QRS-Komplex zu berechnen. Die Gr6IJe dieser Felder spiegelt die Energie 
5 wider die wShrend der Anregung der Vorkammern und Ventriculi des Herzens erzeugt 
wird. Femer kann das VerhSltnis dieser Felder berechnet werden, um so die elektri- 
sche Aktivitat der Vorkammern im Vergleich zu elektrischen AktivitSt der Ventrikel ab- 
zuschatzen. Es hat sich gezeigt, da(i dieser Parameter mit dem Grad des Herzversa- 
gens in Beziehung steht. 

10 Auf derselben Analysestufe kdnnen die Daten eines Vektor-MCG hinzugezogen 

werden, die in speziellen einpositionalen Leads-Systemen registriert wird. Die Amplitu- 
denwerte der Y- und Z-Komponenten wahrend der P-Welle, dem QRS-Komplex 
und dem ST-T-lntervall, die Amplitude und Richtung des raumlich maximalen QRS- 
Vectors, die QRS- Vector-Da uer usw. konnen kalkuliert und fur fur die Diagnose myo- 

15 kardialer Ischamie verwendet werden. Ubiichenweise werden die Analysen zusammen 
mit einer EKG-Analyse durchgefuhrt. 

Der Informationswert der MCG-Daten wSchst dramatisch mit dem Obergang zu 
den folgenden Stufen der Analyse, insbesondere zum Kartographieren des magnetic 
schen Feldes (MF-Mapping). Dieses Verfahren bedeutet die Konstruktion von Vertei- 

20 lungskarten der Induktlon des magnetischen Feldes in MeBpunkten zu bestimmten 
Momenten innerhalb des Kardoizyklus. Diese Karten werden nach dem Prinzip der 
geographlschen Karten konstruiert, d.h. Bereiche mit identischen Werten bestimmter 
Paranrieter besitzen die selbe Farbe. Es ist wichtig zu verstehen, daS jede Karte aus 
alien Me&punkten resultiert. 

25 Die Interpretation der elektrischen Aktivitat des Herzens basierend auf diesem 

Ansatz hat eine Reihe von wesentlichen Vorteilen. ZunSchst kann mit der Hilfe von 
Interpretationsverfahren alle wichtigen Daten. einschliefSlich der Daten. die zwischen 
den Punkten des MeUgitters erhalten werden kSnnten. in die Berechnung einbezogen 
werden. Zweitens stellen Karten eine naturliche Projektion der elektromagnetischen 

30 Phanomene dar, die Qberhalb der Thorax-Oberfiache an verschiedenen Stellen des 
Herzens. wenn auch nur ungefahr, registriert wurden. 

Vorteilhaft kdnnen zwei besonders wichtige Eigenschaften jeder magnetischen 
Feldverteilungskarte automatisch bestimmt werden: erstens die Anzahl der Extrema 
des magnetischen Feldes (im physikalischen Sinn sind lokale Extrema eines magneti- 

35 schen Feldes Punkte mit maximalen Werten im Vergleich zu benachbarten Punkten). 
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mit anderen Worten also die InhomogenitSt der Karte, und zweitens die wechselseitige 
Anordnung dieser Extrema. 

Die Homogenitat der magnetischen Feldkarte spiegelt die Homogenitat der 
elektrischen Quelle, die diese Karte induziert. wieder. Dies wiederum zeigt, dali es 
5 keine myocardischen Orte gibt, die sich wesentlich von benachbarten Zonen hinsicht- 
lich ihrer elektrophysiologischen Eigenschaften unterscheiden, so da& es keine lokalen 
Verietzungsstrome gibt. Normalerweise hat die Karte zu jedem Zeitpunkt des Kardio- 
zyklus eine Dipol-Struktur. d.h., es existiert nur ein Minimum und ein Maximum. Es ist 
kiar, da(i das Auftreten zusatzlicher Extrema das Vorhandensein zusatzlicher lokaler 

10 Str5me belegt. 

Aus physiologischer Sicht ist eine weitere wichtige Charakteristik der Karten. 
wie bereits oben erwahnt, die wechselseitige Anordnung der Minima und Maxima des 
magnetischen Feldes. Falls man eine Linie von einem Minimum zu einem Maximum 
und dann wiederum diese Linie entgegen dem Uhrzeigersinn um 90'* dreht, so wird die 

15 erhaltene Richtung der Orientierung des entsprechenden Ladungsdipols entsprechen. 
Diese Orientierung reflektiert die Richtung der Ausbreitung der Anregungsfront in dem 
jeweils betrachteten Moment des Kardiozyklus. 

Basierend auf diesen Prinzipien werden die Karten von bestimmten charakte- 
ristischen Zeitpunkten des Kardiozyklus. z.B. des QRS-Einsatzes, R-max. QRS-offset. 

20 T-max. T-offset» visuell analysiert. Es kOnnen integrierte Karten, die wahrend des ge- 
samten QRS-Komplexes und/oder dem St-T-lntervall berechnet wurden. untersucht 
werden. Die qualitativen visuellen Analyseverfahren fur magnetische Feldkarten sind 
ausreichend, um eine allgemeine Darstellung wichtiger Eigenschaften des elektrophy- 
sikaiischen Prozesses im Myocardium in jedem individuellen Fall zu erhalten, aber sis 

25 sind nicht in der Lage, eine quantitative Beschreibung der aufgedeckten Eigenschaften 
zu geben und eriauben es nicht, statistische Parameter fur eine Gruppe von Patienten 
zu erhalten. Deshalb ist derfolgende Schritt der Analyse magnetischer Feldvertei- 
lungskarten die Anwendung der quantitativen Kriterien. 

Es wurden verschiedene quantitative Kriterien fur einen vielseitige AbschStzung 

30 verschiedener Karten und Zeitserien von Karten entwickelt. 

Der einfachste Ansatz ist die Berechnung der Anzahl von Extrema in jeder 
Karte und in Ausweitung auf alle untersuchten Intervalle von Kardlozyklen. Dabei wird 
auch der relative "Glattheits-lndex" angewandt, der die Summe der Korelationsfaktoren 
zwischen vier aufeinanderfolgenden Karten am Anfang des ST-Segmentes darstellt. 
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Daneben ist auch das auf der AbschStzung der Komplexitat der Trajektorien der 
Extrema wShrend der ventricularen Anregung basierende Kriterium bekannt. Als weite- 
res quantitatives Kriterium kann die Variabllitat des VerhSltnisses der grdlSten positiven 
zu den grOBten negativen Extremwerten wShrend des ST-T-lntervalls verwendet wer- 
5 den. Ferner ist ein Homogenitats-Koeffizient bekannt, der zur integralen AbschStzung 
der Anzahl von Extrema und deren Scharfe uber das ST-T-lntervall abzuschatzen. Ein 
interessanter Ansatz besteht in einer speziellen r^umlichen Transformation (KLM- 
Transformation) der magnetischen Feldverteilungskarten und der Berechnung der 
nicht-dipolaren Beitrage in jeder Karte. Manchmal warden auch andere quantitative 

10 Parameter venA^endet. Daneben wurden einige andere Ansatze der magnetischen 
Kartographie in verschiedenen AufsStzen beschrieben. So wurde in einem Aufsatz 
simultan mit den traditionellen Karten einer magnetischen Feldverteilung Karten ver- 
wendet, die eine raumliche Verteilung eines quantitativen Indexes darstellt. der die 
Zerlegung eines QRS-Komplexes in jedem Punkt eines MeBgltters widerspiegelt. Dies 

15 erhoht die Genauigkeit der Untersuchung bei Patienten mit VT. Die raumliche Vertei- 
lung des QT-lntervalls und die Dauer seiner Teile verbunden mit der Berechnung eini- 
ger Glattheits-lndizes wird in einem anderen Aufsatz beschrieben. Der Zweck all dieser 
verschiedenen Indizes ist der selbe, nSmlich eine quantitative Inhomogenitatsgradab- 
sch§teung der magnetischen Feldverteilungskarten und damit - bis zu einem gewissen 

20 Grad - der Welle, die dieses Feld produziert, anzugeben. 

Dabei ist es jedenfalls notwendig. zu verstehen. daB die AbschStzung der elek- 
trischen PhSnomene im Myocardium auf Basis der magnetischen Feldverteilung keine 
direkte Abschgtzung ist. 

Eine mehr unmittelbare und physiologisch wertvollere Infomiation kann aus der 

25 MCG-Analyse auf Basis der LOsung des inversen elektrodynamischen Problems ge- 
wonnen werden. Losung des inversen elektrodynamischen Problems betreffend die 
Kardiologie meint die Rekonstruktion der elektrischen Ereignisse im Herzen auf Basis 
der auf der Oberfldche eines menschlichen Kdrpers ausgefQhrten Messungen. Im Fall 
der MCG wird die Messung nicht auf der OberflSche eines Korpers, sondern oberhalb 

30 der Oberfldche in einer Medebene ausgefQhrt. 

Es ist grundsatzlich unmSglich, die Eigenschaften elektrischer Quellen im Her- 
zen. die das EKG erzeugen, allein auf der Basis dieses EKGs zu berechnen. Es wurde 
bereits gesagt, dali die LeitfShigkeit von Geweben und die Form eines Korpers einen 
wesentlich geringeren Einflu(i auf das MCG (verglichen mit dem EKG) besitzt. so dali 

35 die raumliche Auflosung von MCG wesentlich grolier ist als die von EKG. Deshalb er- 
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laubt ein Magnetokardiogramm die wesentlich prazisere LOsung des inversen Pro- 
blems, obwohl die auf Basis eines MCG rekonstruierten Quellen mit Sicherheit ideali- 
siert und nicht komplett mit den tatsachlichen ubereinstimmen. 

Es gibt wenigstens drei Ebenen der Nachbildung und Darstellung bioelektri- 
5 scher Quellen im Herzen. 

Die erste Ebene ist eine Darstellung einer Quelle als aquivalenter Dipot. Dabei 
wird angenommen, dali die gesamte elektrische Aktivitat des Herzens in einem Punkt 
fokusiert ist. Eine solche Darstellung bedeutet nicht, daK das Herz tatsdchlich eine 
Punktquelle ist. Es bedeutet. dad die Resultate seiner AktivitSt auf der OberfiSche ei- 

10 nes K6rpers Equivalent zu den Effekten sind, die gemessen werden kOnnten, wenn 
eine Punktquelle prasent ware. Eine solche Darstellung der Quelle dient als Warneh- 
mungsbasis einer Vector-Kardiographie. Es ist klar. dad es nicht eriaubt, die eigenen 
Aktivitaten verschiedener Telle des Herzens zu bestimmen. 

Die zweite Ebene der Datendarstellung auf Basis der Losung des inversen 

15 Problems ist die Rekonstruktion der Quellen In Form von Ladungsvertellungen in einer 
Schicht. Dabei kdnnen zwei Ansatze venwendet werden. Der erste Ansatz besteht in 
der Interpretation einer magnetischen Feldquelle als eine Karte der Verteilung von La- 
dungsdichtevektoren, der zweite Ansatz erIaubt es, eine Karte von festen Ladungslini- 
en zu zeichnen und ist vielversprechender. Da Bild der Ladungsverteilung eriaubt be- 

20 reits die simultane AbschStzung der Charakteristiken verschiedener Quellen und der 
Anregung verschiedener Telle des Herzens. 

SchlieSlich besteht die dritte Ebene in der Rekonstruktion einer rSumlichen, 
dreidimensionalen bioelektrischen Quelle, d.h. die Wiederherstellung von Quellen die 
der RealitSt am nSchsten kommen. Es ist jedoch klar, dad die Rekonstruktion von 

25 dreidimensionalen Quellen die Anwendung hdchstkomplizierter physikalischer Modelle 
und mathematischer Algorithmen erfordert. 

Die Bestimmung der Position einer Punktquelle im Moment des Anfangs eines 
ectopischen QRS-Komplexes wird angewandt um den Ursprungsort ventrikularer Ar- 
hythmie zu bestimmen. und das selbe Verfahren wird in der Delta-Welle zur zusStzli- 

30 Chen Lokalisation des Pfades angewandt. Die Starke des effektiven Dipols im R- 
Maximum-Punkt wurde auch angewandt als Kriterium zur Abschatzung des Rislkos 
von AbstoHungen nach Herztransplantationen. 

Andere Parameter des effektiven Dipols - seiner Orientierung in einer Fron- 
tebene in charakteristischen Momenten des Kardiozyklus - werden als Kriterium fQr 

35 Patienten bei der Ischemie-Diagnose venA^endet. In unserer Sicht kOnnen sehr wertvol- 
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le Informationen aus der AbschStzung des Dipol-Verhaltens nicht nur in bestimmten 
charakteristischen Momenten des Kardiozyklus. sondern wahrend des gesamten An- 
regungsprozesses gewonnen werden. So ist z.B. die Form der Kurven der Dipol- 
Orientierung und die Starke wShrend des QR-lntervalls ein empfindiicher Marker fur 
5 Ml. In unseren AufsStzen analysieren wir die Kun/en aquivalenter Dipol-Tiefen wSlirend 
der ventrikularen Anregung, d.h. den Charakter der Dipol-Bewegung in einer Vonwarts- 
Ruckwartsrichtung. Wir gehen davon aus, daft zu einem gewissen Grad die Kurve die 
Verteilung der Quellen in einer sagittalen Ebene wiedergibt, was ubiicherweise aufter- 
halb der Analyse verbleibt, da MCG ubiicherweise nur in einer frontalen Ebene regi- 

10 striert wird. Deshalb ist es moglich. die gemeinsame Analyse von Feld-(oder La- 

dungs-)Verteilungskarten und die Kurve von einer Quellentiefe als vereinfachte Analy- 
se einer dreidimensionalen Ladungsverteilung anzusehen. Wie erwartet erscheint die 
Analyse sehr informativ fur die Diagnose von Ml einschlieftlich Nicht-Q Ml zu sein. 

Die Analyse der Verteilungskarten der Ladungs-Dichtevectoren Oder Ladungs- 

15 linien gibt zusatzliche Optionen im Vergleich zu den vorangehenden Analyseebenen. 
Es ist wichtig zu bemerken. dad diese Niveaus der Analyse charakteristisch nur fur 
MCG sind und. sofern wir wissen, keine Analogie in der EKG-Diagnose hat. Sie er- 
moglicht es, ein Bild der Ladungsverteilungen im Myocardium einschlieillich der soge- 
nannten •VerletzungsstrOme" zu visualisieren. 

20 Anwesenheit und Richtung dieser Str6me wird visuell abgeschatzt Als sine 

quantitative Abschatzung der Stromverteilung wird der Stromdichtewert venvendet, 
z.B. in der Spitze des QRS-Komplexes. Wir analysieren die Kurven der maximaien und 
der durchschnittlichen Stromdichte in jeder Karte wahrend des gesamten Kariozyklus. 
Eine solche Analyse erscheint sehr brauchbar fOr die verschiedenen Formen der 

25 Ischemie. 

ErfindungsgemaiS wird femer ein automatsiches Verfahren vorgeschlagen . dafi 
dem Arzt die Moglichkeit zur Verfugung stellt, die Anzahl, Richtung, Intensitat und Gr6- 
Be von Stromwirbein In jedem Moment des Kardiozyklus abzuschatzen und damit ihr 
Verhalten wahrend der Depolarisation Oder Repolarisation. 

30 Zweckmaftig ist es, parallel zur Magnetokardiographie es notwendig ist, ein 

Herz-Visualisierungsverfahren anzuwenden, wie z.B. Rontgen, Computertomographie, 
Kemspinresinanztomographie usw,. Welches dieser Verfahren gewahit werden sollte 
hangt von den Anforderungen an die detaillierte Auflosung der anatomischen Informa- 
tionen ab, die wiederum von der konkreten klinischen Aufgabe abhangt. Die Erfahaing 

35 zeigt, daft es in den meisten Fallen ausreicht, einfaches Rontgen anzuwenden Oder 
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daH es sogar mflglich ist, uberhaupt keine Herz-Visualisierungsverfahren anzuwenden, 
ohne dad es zu einer Reduzierung des Wertes der MCG-Untersuchung kame. 

Bei der Analyse gemaii Tabelle 1 ist jede aufeinanderfolgende Stufe der Analy- 
se spezifiziert und entwickelt weiter die Informationen, in der vorangehenden Stufe 
5 erhalten wurden. Beispielsweise wSre die Sensitivitat des Algorithmus bei der Diagno- 
se von IHD, wenn man nur die Analyse der magnetischen Feldverteilungskarten an- 
wendete, 86 %. Falls man zusatzlich dazu die Analyse basierend auf der L6sung des 
inversen Problems anwendete, erhdhte sich die Sensitivitat auf 94 %. 

Die Form der medizinischen Schlulifolgerung, die von dem Magnetokardiogra- 

10 phologen gemacht und an den Kiiniker weitergeleitet wird, ist sehr bedeutend. Es wer- 
den Schlulifolgerungen verwendet, die aus zwei Teilen bestehen. Im ersten Teil wer- 
den quantitative und halb-quantitative Charakteristiken des aktuellen MCG angegeben 
(Homogenitat der Karten, Richtung des ECD und der Vectoren der Stromdichte, 
Stromdichtewerte usw.). 

15 Im zv^eiten Teil der Schlulifolgerung werden die entdeckten Veranderungen mit 

physiologischen Schlussen Oder Typen von Herzkrankheiten automatesch in Bezie- 
hung gesetzt, z.B.: "... Storungen, die wahrend der ventrikularen Repolarisation wur- 
den, zeugen von hoher (mittlerer, geringer) Wahrscheinlichkeit von IHD" Oder"... von 
einem hohen Risiko des Auftretens von Arlthmlen" oder "im Vergleich mit vorangehen- 

20 den MCG zeigen sich signifikante positive Veranderungen, die die Effizienz einer anti- 
anginalen (antiarithmischen) Therapie bestatigen". 

Sicher ist das vorgeschlagene Interpretationsschema und die Form der SchlQs- 
se nicht komplett oder das einzig mdgliche. Jedes Interpretationsniveau und die sepa- 
raten Parameter erfordern eine detailtierte Analyse aus Sicht der Informationstheorie. 

25 Nach allem wird es mdglich. einen Satz von formalisierten Entscheidungsregein fQr die 
Diagnose verschiedener Herzzustande auf Basis MCG-Daten anzugeben. 

Es steht ohne Zweifel, dafi die Entwicklung und Standardisierung eines Sche- 
mas fUr die medizinische Interpretation von MCG eine der bedeutendsten notwendigen 
Bedingungen fur eine breite Anwendbarkeit des Verfahrens ist. 

30 Jedenfails ist es bereits jetzt moglich, anzuerkennen, dall die geschaffenen 

Kriterien und die Interpretationsschemata die Anwendung von MCG als reales und 
wirklich nutzliches diagnostisches Werkzeug anzuerkennen. 

Im allgemeinen k6nnen alle Kriterien in vier Gruppen eingeteilt werden. Inner- 
haib der Gruppen sind sie in der Reihenfolge aufsteigender Analysestufen angeordnet. 

35 Gruppe 1 . Kriterien zur Abschatzung des Signal-Rausch-Verhaitnisses. 
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A. Visuelle Abschatzung der Hochfrequenz-Niedrig-Amplitutde-Wellen entlang 
durchschnittlicher MCG-Kurven. 

Gruppe 2. Kriterien fur die Abschatzung der Anregungshomogenitat. 

A. Visuelle Abschatzung der Verschiedenheit der magnetischen Feldkarten. 
5 B. Quantitative Analyse der Anzahl von Extrema. Dabei ist der Normalwert - 

nicht mehr als drei Extrema in jeder Karte. 

C. Koeffizient der Homogenitat (CH). Der Normalwert ist nicht mehr als 0,95. 
D*. 5mi„/nax wShrend des ST-T-Intervalls (min/max-Vcrhaitnis der maximalen negativen 
Extremwerte zu den minimalen negativen Extremwerten zu jedem Zeltpunkt). Der 
10 . Normalwert ist nicht grQIier als 0,15. 

E*. 6Zo wahrend ST-T-lntervall (Zo - Tiefe der effektiven Dipole in jedem Zelt- 
punkt). Der 

Normalwert ist nicht mehr als 0,20. 

F. Visuelle Abschatzung der Ungleichheit der Stromverteilungskarten. 
15 G*. Quantitative Analyse der Stromwirbelanzahl. Der Normalwert ist nicht mehr 

als drei Wirbel in jeder Karte. 

Die hauptsachliche klinische Signifikanz der Gruppe, von Kriterien ist die Be- 
stimmung von Herzstorungen im allgemeinen und im besonderen die Bestimmung des 
Risikos von kardialer Arithmie und die Abschatzung der Effektivitat einer Antiarithmi- 
20 schen Therapie. 

Gruppe 3. Kriterien fQr die Analyse der Anregung der Richtung der Wellenfron- 
tausbreitung. 

A. Visuelle Abschatzung der ungefahren Richtung der entsprechenden Strom- 
dipole basierend auf den magnetischen Feldverteilungskarten. 

25 B*. Quantitative Analyse der entsprechenden Stromdipolrichtung basierend auf 

den magnetischen Feldverteilungskarten. Die Normale Richtung wahrend des ST-T- 
Intervalls ist nach links unten (wahrend nicht mehr als 2/3 der Dauer des ST-T- 
Intervalls). Die Nomiale Richtung wahrend des QRS-Komplexes besteht aus drei Pha- 
sen: 1. Phase nach rechts unten, 2. Phase nach links unten, 3. Phase nach oben. 

30 C*. Analyse der effektiven Dipoltiefe. 

Normaler Parameter wahrend des ST-T-lntefvalls - die Tiefe im J-Punkt sollte haupt- 
sachlich in dem Interval! tiegen. 

Normaler Parameter wahrend des QRS-Komplexes - vier unterschiedliche Bewegun- 
gen eines Dipols in der Tiefe existieren. Die erste Bewegung ist vorv^arts gerichtet, die 
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zweite Bewegung (die haupts^chliche) ist rOckwarts gerichtet, die dritte Bewegung ist 
vorwSrts gerichtet, die vierte Bewegung ist ruckwarts gerichtet. 

D. Visuelle Abschatzung der Richtung der Stromausbreitung basierend auf den 
Karten der Stromlinien und den Karten der Stromdichtevectoren. 
5 E*. Quantitative Analyse der Ausbreitungsrichtung des Stroms basierend auf 

den Karten der Stromlinien und den Karten der Stromdichtevectoren. 
Normale Richtung w^hrend des ST-T-lntervalls ist nach links unten (wahrend nicht 
mehr als 2/3 der Dauer des ST-T-lntervalls). 

Normale Richtung wahrend des QRS-Komplexes besteht aus drei Phasen: 1. Phase 
10 ' nach rechts unten gerichtet, 2, Phase nach links unten gerichtet, 3. Phase nach oben 
gerichtet. 

Die hauptsachliche klinische Signifikanz dieser Gruppe von Kriterien ist die 
Bestimmung von Herzstorungen im allgemeinen. insbesondere die Diagnose verschie- 
dener Formen der Ischemie und die Abschatzung der Effizienz einer antiischemischen 
15 Therapie. 

Gruppe 4. Kriterien fQr die Analyse absoluter und relativer Intensitat des elektri- 
schen Prozesses im Myocardium. 

A*. Analyse der relativen Vorkammern zur ventrikularen elektrischen Aktivitat 
(P/QRS-lntegralverhaltnis). Der 
20 Normalwert - nicht mehr als 0, 1 3. 

B*. Analyse der Stromdichtewerte. 
Normalwerte sind: 

- das Verhaitnis der Stromdichte im Zeitpunkt 80 msec nach ST-Einsatz zur Stromdich- 
te im J-Punkt sollte nicht kleiner als 2.5 sein 

25 - das Verhaitnis der Stromdichtewerte am R-max-Punkt zu dem am T-max-Punkt sollte 
nicht grd&er als 3,5 sein. 

Die hauptsachliche klinische Signifikanz dieser Kriteriengruppe ist die Bestim- 
mung von Herzstorungen im allgemeinen, insbesondere die Diagnose verschiedener 
Formen der Ischemie und des Herzversagens, sowie die Abschatzung der Therapieef- 

30 fizienz. 

Die groUe Anzahl und Struktur der Kriterien spiegelt die vielen Aspekte der im 
MCG enthaltenen Informationen und auch die historische Richtung der Softwareent- 
wicklung wieder: von der morphologischen Analyse der MCG-Kurven zur Lfisung des 
zweidimensionalen inversen Problems. 
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Jetzt werden in der praktischen Arbeit alle oben genannten Satze verwendet. 
Manchmal gibt es bei der Analyse des MCG von einigen Patienten Inconsistenzen zwi- 
schen diesen Kriterien (beispielswelse entspricht ein Indilcator der Norm, wShrend ein 
anderer einen abnormalen Wert zeigt). Das ist der Grund, warum es notwendig ist. 
5 basierend auf disl<renriinater Clusteranalyse Oder anderen mathematischen Methoden 
fur die Mustererkennung sorgfaltig die informativen MSglichkeiten jedes Indikators zu 
bestimmen und sodann einen Satz von fomnalisierten Losungsregein fQr verschiedene 
klinische Aufgaben anzugeben, die automatisch arbeiten kdnnen. Dann wird es mdg- 
lich, einige Kriterien zuruckzuweisen. 
10 Diagnostische Kriterien fur die Diagnose myocardiaier Ischemie unter Verwen- 

dung magnetocardiographischer Daten. 
Qualitative Analyse der Karten der Verteilung des ventrikularen magetischen 
Repolarisationsfeldes 

gemali der Rechte-Hand-Regel wird die ungefahre Richtung der entsprechen- 
15 den Stromdipole in jeder Karte bestimmt und ihr Abweichen von der normalen Richtung 
nach links unten ermittelt. 

die Anwesenheit zusatzlicher Extreme in jeder magnetischen Feldverteilungs- 
karte wird ermittelt. 

Quantitative Analyse der ventrikularen magnetischen Repolarisationsfeldvertei- 
20 lungskarten. 

Die Ordnungszahl der Karte wird bestimmt ausgehend von derjenigen. in wel- 
cher die Richtung des Squivalenten ventrikularen Repolarisations-Stromdipols stabil 
normal wurde, d.h. nach links unten zeigte (wahrend nicht mehr 1/3 der ST-T- 
Intervalldauer, wobel je grd&er die Ordnungszahl der genannten Karte desto gr5&er die 
25 Schwere der Ischemie). 

- * Dauer der Existenz zusatzlicher Extreme wdhrend des Repolarisationsprozes- 
ses wird abgeschatzt (nicht mehr als 1/3 der Dauer des ST-T-lntervalls. je grdBer die 
Dauer der Existenz zusatzlicher Extreme, desto groiier die Schwere der Ischemie). 
Analyse der effektiven Dipol-Parameter innerhalb des ST-T-lntervalls 
30 - Die effektive Dipol-Tiefe in J-Punkt wird abgeschatzt (die Tiefe in diesem Punkt 
sollte die grolite innerhalb des ST-T-lntervalls sein). 

Qualitative Analyse der ventrikularen Repolarisations-Spannungslinien und 
Dichtevectorkarten. 

Die ungefahre Richtung der meisten stromdichtevectoren in jeder Karte wird 
35 abgeschatzt und ihr Abweichen von der normalen Richtung nach links unten bestimmt. 
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Die Anwesenheit zusatziicher Stromwirbel in jeder Stromlinienkarte wird ermit- 

telt. 

Quantitative Analyse der vetrikularen Repolarisations-Spannungslinien und 
Dichte-Vectorkarten 

5 - Die Ordnungszahl der Karte wird bestimmt, ausgehend von der Karte, in wel- 
cher die Richtung der meisten stromdicfitevectoren stabil normal, d.h. nach links unten 
gerichtet, wurde (nicht mehr als 1/3 der Dauer des ST-T-lntervalls, je grolJer die Kar- 
tenordnungszahl, desto gr6(Jer die Schwere der Ischemie). 

Die Dauer der Existenz zusatziicher Spannungswirbel wird ermittelt (nicht mehr 
10 als 1/3 der Dauer des ST-T-lntervalls, je grolier die Dauer der Existenz zusatziicher 
Spannungswirbel, desto grfiBer die Schwere der Ischemie). 

Quantitative Stromdichte Parameter innerhalb des ST-T-intervails 
Das Verhaltnis der Stromdichte im Zeitpunkt 80msec nach ST-Einsatz zur 
Stromdichte im J-Punkt wird bestimmt (es sollte nicht kleiner als 2,5 sein). 
1 5 - Das Verhaltnis der Stromdichtewerte im R-max-Punkt zu dem im T-max-Punkt 
wird berechnet (es sollte nicht grSlier sein als 3,5). 

Jede nachfoigende Stufe der Analyse erganzt die Informationen aus der voran- 
gehenden Stufe und weitet diese aus. 

Auf der Basis aller oben genannten Kriterien wird der SchluH Uber das Vorlie- 
20 gen und die Schwere der myocardialen Ischemie gezogen. 

Es werden folgende vorteilhafte Weiterbildungen vorgeschlagen. Zwei von Ih- 
nen sind die grundlegenden und offenbaren einen "Instrumentensatz" fur die umfang- 
reiche Analyse der ventrikularen Repolarisatlon bzw. Depotarisation. Zwei weitere Pa- 
tente sind der Unterscheidung der bedeutendsten Herzkrankheiten gewidmet: chronl- 
25 sche Ischemie und dem myocardialen Infarkt. In der Zukunft wird es moglich sein, dia- 
gnostische Verfahren fur einige andere Herzkrankheiten anzugeben. 

1 ) Verfahren zur Evaluation der ventrikularen Repolarisatlon unter Venvendung 
magnetocardiographischer Daten 

3) Das Problem ist das folgende: die Bestimmung der elektrophysikalischen 
30 Eigenschaften der ventrikularen Repolarisatlon, die die grolite Signifikanz besitzen, 
unter Verwendung der Magnetocardiographie, wobei die Eigenschaften als Kriterium 
fur die Separation von normalen und pathologischen Funktionszustanden des Herzens 
und zur Formulierung von Schlussen uber bestimmte Herzkrankheiten dienen konnen. 
Das vorgeschlagene Verfahren basiert auf verschiedenen aufeinanderfolgen- 
35 den Stufen der MCG-Analyse: visuelle qualitative und quantitative Analyse der magne- 
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tischen Feldverteilungskarten. Analyse der effektiven Stromdipollokalisation. qualitative 
und quantitative Kriterien der Stromverteilung. All diese Schritte eriauben es, eine um- 
fassende, akkurate und vielseitige Abschatzung der Homogenitat der Anregung, der 
Richtung der Strome, der Leistungscharaktertistik der Anregung in jedem Zeitpunkt und 
5 das Verhalten all dieser Parameter wahrend des gesamten Repolarisationsprozesses 
auszufuhren. Neben diesen Kriterien eriaubt es die Analyse der effektiven Dipoltiefe 
nicht nur zweidimensionale Verteilungen einer Quelle, sondern. bis zu einem gewissen 
Grad, ihre dreidimensionale Verteiiung zu erhalten. 

Wesentlich neue Innovationen der angebotenen Erfindung sind: 
10 a) aufeinanderfolgende Anwendungen verschiedener Stufen der MCG-Analyse 

b) Abschatzung der Richtungsanderungen des aquivalenten Stromdipols auf den ma- 
gnetischen Feldverteilungskarten und auf den Karten der Stron^linienverteilung wah- 
rend der Repolarisation. 

c) Abschatzung der effektiven Dipoltiefenkun/e. 

15 d) Auf einer Abschatzung der Leistungscharakteristik basierende Kriterien (Kurven 
maximaler und durchschnittlicher Dichte eines Stroms wahrend der Repolarisation). 

Alle oben genannten Kriterien eriauben es, den ventrikularen Repolarisations- 
prozeH umfassender, detaillierter und genauer abzuschatzen und auf dieser Basis 
SchlQsse uber den funktionalen Zustand des Myocardiums zu Ziehen. 

20 Verfahren zur Evaluation der ventrikularen Depolarisation unter VenA/endung 

magnetocardiaographischer Daten. 

Das Problem ist das folgende: Bestimmung der signifikantesten elektrophysika- 
lischen Eigenschaften der ventrikularen Depolarisation unter Verwendung der Magne- 
tocardiographie, die als Kriterien fur die Unterscheidung von nomnalen und Pathologi- 

25 schen Funktionszustanden des Herzens und auch zur Gewinnung von Aussagen Uber 
verschiedene Herzkrankheiten dienen k6nnen. 

Es gibt verschiedene Verfahren zur Evaluation der ventrikularen Depolarisation, 
namlich: Verfahren der spaten Felder. Verfahren der visuellen Abschatzung momenta- 
ner magnetischer Feldverteilungskarten und zeitlich integrierter Karten, Verfahren der 

30 qualitativen und quantitativen Evaluation von Residuen-Karten, Verfahren der Evalua- 
tion der effektiven Stromdipol-Parameter, Verfahren der Stromdichteberechnung. 

Das erste dieser Verfahren reprasentiert die Spektralanalyse des letzten Teils 
des QRS-Komplexes, vergleichbar mit der Analyse der spaten Potentiale beim EKG. 
Das zweite Verfahren stellt die Analyse der wechselseitigen Anordnung von positiven 

35 und negativen Extrema in jeder magnetischen Feldverteilungskarte. Das dritte Verfah- 
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ren stellt einer visuelle und quantitative Abschatzung der Unterschiede zwischen der 
aktuellen Karte und einer rekonstruierten normalen Karte dar. Das vierte Verfahren 
stellt eine Abschatzung der Orientierung und Starke der Stromdipole wShrend des QR- 
Intervalls dar Das fOnfte Verfahren stellt eine Abschatzung der Stromdichteverteilung 
5 und der Werte in den QRS-maximum-Vectoren zu den verschiedenen Zeitpunkten der 
Depolarisation dar. 

Jedes dieser bekannten Verfahren eriaubt die Abschatzung nur eines Aspektes 
der ventrikularen Depolarisation: Homogenitat der Anregung oder Richtung der aquiva- 
lenten Stromdipole oder Stromdichte in definierten Zeitpunkten der Depolarisation. 

10 Deshalb gibt es keine quantitativen Kriterien, die es eriaubten, den Depolarisationspro- 
zefi nicht nur in diskreten Zeitpunkten, sondern wahrend des gesamten Prozesses 
abzuschatzen. Jedes dieser Verfahren basiert nur auf einer Stufe der Analyse: auf den 
magnetischen Feldverteiiungskarten oder den Parametern effektiver Dipole oder den 
Stromdichte- Verteilungskarten. 

15 Das vorgeschlagene Verfahren basiert auf verschiedenen aufeinanderfolgen- 

den Stufen der MCG-Analyse: visuelle qualitative und quantitative Analyse der magne- 
tischen Feldverteiiungskarten, Analyse der Lokalisatlon effektiver Stromdipole, qualita- 
tive und quantitative Kriterien der Stromverteilung. All diese Schritte eriauben es. eine 
umfassende, genaue und vielseitige Abschatzung der Homogenitat der Anregung, der 

20 Richtung der Strbme, der Leistungscharakteristik der Anregung in jedem Zeitpunkt und 
das Verhalten all dieser Parameter wahrend des gesamten Depolarisationsprozesses. 
Daneben eriauben es die auf der Analyse der effektiven Dipoltiefe basierenden Kriteri- 
en, nicht nur die zweidimensionalen Verteilungen einer Quelle, sondern, bis zu einem 
gewissen Grade, auch deren dreidimensionale Verteilungen zu erhalten. 

25 Wesentliche neue Innovationen der angebotenen Erfindung sind: 

a) konsekutive Anwendung verschledener Stufen der MCG-Analyse. 

b) Abschatzung der Richtungsanderungen der aquivalenten Stromdipole auf magneti- 
schen Feldverteiiungskarten und auf Karten der Stromlinienverteilung wahrend der 
Depolarisation. 

30 c) Abschatzung der effektiven Dipoltiefenkurven. 

d) Kriterien, die auf einer Leistungscharakteristikabschatzung (Kurven der maximalen 
und der durchschnittlichen Dichte eines Stroms wahrend der Depolarisation) wahrend 
des gesamten Depolarisationsprozesses beruhen. 
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Alle oben genannten Kriterien erlauben es, den ventrikularen Depolarisations- 
prozeR umfassender. detaillierter und genauer abzuschatzen und auf dieser Basis 
Schlusse iiber den funktionalen myocardialen Zustand zu Ziehen. 

Verfahren zur Diagnose myocardialer Ischemie unter Verwendung magneto- 
5 cardiographischer Daten be! Patienten mit nicht-verandertem Ruhe-EKG. 

Das Problem is! das folgende: MCG-Bestimmung der elektrophysikalischen 
Eigenschaften der ventrikularen Repolarisation, die von myocardialer Ischemie her- 
rOhrt. Dies ist besonders bedeutend bei Patienten mit nicht-infonnativen EKG. 

Es gibt zwei Verfahren zur MCG-lschemie-Dlagnose bei Patienten mit nicht- 

10 verandertem EKG, namlich: Verfahren der visuellen AbschStzung magnetischer Feld- 
verteilungskarten und Evaluationsverfahren der effektiven Stromdipolparameter. Das 
erste dieser Verfahren stellt eine visuelle Abschatzung der Quantitat und Anordnung 
der magnetischen Extrema und der Richtung der Dipole in verschiedenen Momenten 
wahrend der Repolarisationsperiode dar. 

15 Das zweite Verfahren stellt eine Abschatzung der Richtung und der Starke der 

Stromdipole In definierten Momenten des Kardiozyklus dar. 

Nachteil der ersten dieser Verfahren ist die Tatsache. dalS nur die rSumliche 
Verteilung des magnetischen Feldes erhalten werden kann ohne BerQcksichtigung der 
Richtung der Dipole und der Leistungsckarakteristik. Umgekehrt berucksichtigt das 

20 zweite Verfahren nicht die rSumliche Verteilung eines magnetischen Feldes. Gemein- 
samer Nachteil dieser beiden Verfahren ist das Fehlen quantitativer Kriterien, die es 
eriaubten, den RepolarisationsprozeB nicht nur in einem diskreten Zeitpunkt, sondern 
wShrend des gesamten Prozesses abzuschatzen. Diese Verfahren sind nicht in der 
Lage, die Schwere der Ischemie abzuschatzen. Jedes von ihnen eriaubt es lediglich, 

25 nur einen Aspekt der ventrikularen Repolarisation abzuschatzen: die Homogenitat der 
Anregung Oder die Richtung der aquivalenten Stromdipole und die Starke in bestimm- 
ten Momenten der Repolarisation. 

Das angebotene Verfahren basiert auf verschiedenen aufelnanderfolgenden 
Stufen der MCG-Analyse: visuelle qualitative und quantitative Analyse der magneti- 

30 schen Feldverteilungskarten, Analyse der Lokalisation der effektiven Stromdipole, 

qualitative und quantitative Kriterien der Stromverteilung. Alle diese Schritte eriauben 
es, eine umfassende, genaue und vielseitige Abschatzung der Anregungshomogenitat, 
der Richtung der Stroma, der Leistungscharakteristik der Anregung in jedem Zeitpunkt 
und des Verhaltens all dieser Parameter wahrend des gesamten Repolarisationspro- 

35 zesses. Daneben eriauben es die auf der Analyse der effektiven Dipottiefe basierenden 
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Kriterien nicht nur zweidimensionalen Verteiiungen einer Quelle, sondern, bis zu einem 
gewissen Grad, auch deren dreidimensionale Verteiiungen zu erhalten. Dies macht es 
mOglich. die Sensitivitat und Spezifitat eines Verfahrens signifikant zu verbessern und 
nicht nur das Vorliegen einer Ischemie zu diagnostizieren, sondern auch ihren Grad 
5 abzuschStzen. 

Wesentliche neue Innovationen der angebotenen Erfindung sind: 

a) aufeinanderfolgende Anwendung verschiedener Stufen der MCG-Analyse 

b) Abschatzung der Richtungsanderungen der Squivalenten Stronndipole auf magneti- 
schen Feldverteilungskarten und auf Karten der Verteilung der Stromlinien wahrend 

10 der Repolarisation 

c) Evaluation zusatzlicher Extreme und der Strombereichs-Lebenszeit wShrend der 
Repolarisation. 

d) Abschatzung der effektiven Dipoltiefenkurven. 

e) Kriterien basierend auf einer Abschatzung der Leistungscharakteristiken (Kurven der 
15 maximalen und der durchschnittlichen Dichte eines Stroms wahrend der Repolarisati- 
on). 

Alle oben genannten Kriterien eriauben es, die Sensitivitat und Genauigkeit der 
Diagnose myocardialer Ischennie mittels MCG be! Patienten mit nicht-verandertem 
EKG signifikant zu verbessern. 
20 Verfahren zur Diagnose myocardialer Nekrose unter Verwendung magnetocar- 

diographischer Daten. 

Das Problem ist das folgende: Bestimmung der myocardlalen Nekrose, die von 
einem myocardlalen lnfart<t hen^uhrt. d.h. Entwicklung von MCG-Aqulvalenten der myo- 
cardlalen Nekrose. 

25 Es gibt verschledene Verfahren der Evaluation myocardialer Nekrose unter 

Verwendung von MCG, namlich: Verfahren der morphologischen Analyse des QRS- 
Komplexes, Verfahren der visuellen Abschatzung verzogemngsfreier magnetischer 
Feldverteilungskarten und zeitlich integrierter Karten, Verfahren der qualitativen und 
quantitativen Evaluation von Residuenkarten, Evaluationsverfahren der effektiven 

30 Stromdipolparameter, Verfahren der Berechnung von Stromdichten. 

Das erste Verfahren stellt die Abschatzung des QRS-Komplex-Typus in ver- 
schiedenen MeBpunkten dar, aquivalent zur Standard EKG-Analyse. 

Das zweite Verfahren stellt die Analyse der wechselseitigen Anordnung von 
positiven und negativen Extreme in jeder magnetischen Feldverteilungskarte dar. Das 

35 dritte Verfahren stellt eine visuelle und quantitative Abschatzung der Unterschiede zwi- 
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schen der aktuellen Karte und einer rekonstruierten normalen Karte dar. Das vierte 
Verfahren stellt eine AbschStzung der Richtung und StSrke der Stromdipole wahrend 
des QRS-lnten^alls dar. Das funfte Verfahren stellt eine Abschatzung der Stromdichte- 
verteilung und der Werte in den QRS-Maximum-Vectoren zu verschiedenen Zeitpunk- 
5 ten der Depolarisation dar. 

Jedes dieser Verfahren eriaubt es, nur einen Aspekt der ventrikularen Depolari- 
sation abzuschMtzen: HomogenitSt der Anregung Oder Richtung der Squivalenten 
Stromdipole Oder Stromdichte in definierten Zeitpunkten der Depolarisation. Deshalb 
gibt es keine quantitativen Kriterien. die es eriaubten, den Depoiarisationsprozed nicht 

10 nur in den bestimmten Zeitpunkten, sondern wahrend des gesamten Prozesses abzu- 
schStzen. Jedes dieser Verfahren basiert nur auf einer Stufe der Analyse: auf magnetic 
schen Feldverteilungskarten oder effektiven Dipolparametern oder Stromdichtevertei- 
lungskarten. Nur das zweite Verfahren wird zur Diagnose des Nicht-Q-lnfarktes. also 
bei Patienten mit nicht-informativem oder zweifelhaftem EKG verwandt. 

15 Das vorgeschlagene Verfahren basiert auf verschiedenen aufeinanderfolgen- 

den Stufen der MCG-Analyse: visuelle qualitative und quantitative Analyse der magne- 
tischen Feldverteilungskarten, Analyse der Lokalisation der effektiven Stromdipole, 
qualitative und quantitative Kriterien fQr die Stromverteilung. All diese Schritte eriauben 
eine umfassende. genaue und vielseitige Abschatzung der Homogenitat der Anregung. 

20 der Richtung der StrOme, der Leistungscharakteristik der Anregung in jedem Zeitpunkt 
das Verhalten all dieser Parameter wdhrend des gesamten Depolarisationsprozesses. 
Daneben eriauben es die auf der Analyse der effektiven Dipoltiefe basierenden Kriteri- 
en nicht nur, die zweidimensionale Verteilung einer Quelle, sondern, bis zu einem ge- 
wissen Grade, auch deren dreidimensionale Verteilung zu erhalten. 

25 Die wesentlichen neuen Innovationen der angebotenen Erfindung sind: 

a) aufeinanderfolgende Anwendung verschiedener Stufen der MCG-Analyse. 

b) Abschatzung der RichtungsSnderungen der Squivalenten Stromdipole auf magneti- 
schen Feldverteilungskarten und auf Karten der Stromlinienverteilung wahrend der 
Depolarisation. 

30 c) Abschatzung der effektiven Dipoltiefenkurven. 

d) Kriterien basierend auf der Abschatzung der Leistungscharakteristik (Kurven der 
maximalen und durchschnittlichen Dichte eines Stroms wahrend der Depolarisation) 
wahrend des gesamten Depolarisationsprozesses. 

Alle oben genannten Kriterien eriauben es, die Empfindlichkeit und Genauigkeit 

35 der Diagnose myocardialer Nekrose mittels MCG signifikant zu verbessem. 
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In der Figur 1 sind nach Art eines Flulidiagramm die wichtigstens Schritte bei 
der Durchfuhrung des Verfahrens gemaU Anspruch 1 gezeigt, 

Zur Berechnung der Quellen eines erfaliten biomagnetischen Feldes warden 
zwei parallels z.B. quadratische FiSchen betrachtet, die urn den Abstand h voneinan- 
5 der entfernt sind. Die "obere" FISche wird als MeBebene bezeichnet. die "untere" 
Fiache Ss als Quellenebene. Das magnetische Feld wird durch Gradiometer zweiter 
Ordnung erfalit, dafi uberhalb der MeBebene so plaziert siind, dafi die aufnehmende 
Spule (Pick-up-Coil) in der Mefiebene liegt. Die Messung eines Signals erfolgt an den 
uberiappenden Zentren eines Pick-up-Coils mit n x n Knoten in dieser Ebene. Das 
10 Gitter von m x m Knoten wird auf Sg gewahlt. 

Die Integralgleichung (IE) Ax=y wird gelost, wobei die Koeffizienten einer Matrix 
ay. i=1 , n^, j=1, definiert werden durch den Ausdruck s=Si-2s2+S3. Die Komponenten 
Sr2s2+S3 sind vertikale (z) Komponenten der magnetischen Induktion in den Punkten 
Qi, Q'j. Q"j. wobei Q| ein Knoten des Gitters und die Punkte QV Q" oberhalb von in 
15 den beiden Distanzen b und 2b sind, wobei b die Baselinedes Gradiometers ist. Der 
Koeffizient ay hSngt auch von der Position eines Knotens Mj auf einem Gitter Ss ab. 
Der rechte Tell der Gleichung stellt die Verteilung in den Knoten Qj eines Signals eines 
Gradiometers dar. 

Nach LOsung von IE mit Hilfe der SV-Dekomposition einer Stromfunktion Oder 
20 einer Verteilung der magnetischen Momente ertialt man 5 in einer Quellenebene. Sei- 
ne Isolinien sind Stromlinien. Die Verteilung der Stromdichtevectoren findet man mittels 
der Fonnel j=[grad x. n]. 

Fur die Definition der Abstande zwischen Ebene Ss und S^ - d.h. die 'Tiefe" - 
geht man wie folgt vor: Man lOst das Problem fur die Werte hn, k=1,p und bestimmt 
25 eine Position auf einer Ebene Ss, fOr welche die Norm || x^*\\/ (xOmax minimal ist. Der 

gefundene Wert h k* wird verbessert durch FQhren einer Parabel durch die Punkte h kM. 
h Ks h und Berechnen einer Position eines Minimums dieser Parabel. 
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Stufe der 
Untersuchung 


untersuchtes 

physlkalisches 

Substrat 


informative Kriterien 


Erscheinung der 
resuttierenden information. 


visuelle qualitative 
Untersuchung von 
MCG in jedem von 
36 

Aufnahmepunkten 


nnagnetisches Feld 
gegenOber der 
Zeitabhdnglgkeit 
wShrend des 
Kardiozyklus in jedem 
von 36 

Aufnahmepunkten. 


Signal-Rausch-Verhaitnis. 


Schlusse hinsichtlich der 
Geeignethett von MCG fiir 
weitere Untersuchungen. 


Untersuchung des 
integralen MCG. 


durchschnittliches 
magnetisches Feld 
gegenOber der 
Zettkurve wdhrend des 
Kardiozyklus fUr 
verschiedene 
Aufnahmepunkte. 


Verhditnis der Integrale 
P/QRS. 


Schlusse bezuglich der relativen 
eiektrischen Atria zur 
ventrikularen Aktivltdt. 


visuelle 

(qualitative) 

Untersuchung der 

magnetischen 

Feidvertellungskart 

en. 


unmittelbare und 
rdumllche Karten der 
Verteiiung der 
magnetischen 
Feldintensitdt. 


Auftreten von zusStzlichen 
Orten positiver und 
negativer Werte und deren 
reziproker Position. 


Schlusse bezOglich der 
normalen Oder reduzierten 
Homogenitat des 
elektrophysiologischen 
Prozesses. provisorischer Biick 
in Richtung aufEMF 
(elektromagnetische Kraft - 
electromagnetic force). 


quantitative 
Analyse der 
n^agnetischen 
Feidvertellungskart 
en. 


unmittelbare rSumliche 
Karten der Verteiiung 
der magnetischen 
Fetdintensitdt und 
dessen Variationen 
wahrend des 
Kardiozyklus 


Anzahl der Extrema, IH KH, 
d min/max, Kartennummer 
mit normaler Rtchtung von 
EMF (elektromagnetischer 
Feldstarkevector). 


Schlusse hinsichtlich der 
Homogenitcitsrate der Atria- und 
der ventrikularen Anregung und 
Abnormalitaten der EMF- 
Richtung wShrend der 
ventrikularen Repolarisation. 


Analyse der 
effektiven Punkte 
des Dipolpfades. 


Vahationen der 
effektiven 

Dipolposition wShrend 
des Kardiozyklus. 


Typ der effektiven 
Dipoltiefenkurve wdhrend 
des Kardiozyklus. visuelle 
Kriterien der effektiven 
Dipoipositionsvariabilitat dz 
wdhrend der ventrikularen 
Repolarisation. 


Ausrichtung der Inhomogenitdt 
und anderer 
elektrophysiologischer 
Stdrungsraten und Typus, 
einschlielilich Verteiiung der 
Quellen in der sagittalen Ebene. 


qualitative und 
quantitative 
Analyse der 
zweidimensionalen 
Ladungsverteiiung. 


unmittelbare rdumtlche 
Karten der Stromlinien 
und der Verteiiung der 
Stromdichtevectoren 


visuelle und quantitative 
Kriterien^ die es ertauben. 
die Anzahl der Quetlen, 
deren Ort, relative Intensitat, 
die Richtung der fliefienden 
StrSme und die Variation 
dieser Parameter wShrend 
des Kardiozyklus 
abzuschatzen, Kurven der 
maximalen und der 
durchschnittlichen 
Stromdichte. 


weitere Ausrichtung der 
Inhomogenitdt und anderer 
elektrophysiologischer 
Stcirungsralen und Typus, 
Abschatzung der relativen 
Intensitat der nomnalen und der 
pathologischen Quellen und 
ungefahre Lokalisation der 
ausgefailenen Zonen des 
Myocardiums, Abschdtzung der 
energetischen 
Anregungsparameter 
(Stromdichte). 



integrierte SchlQsse hinsichtlich der Inhomogenrtatsrate und des Typus des myocardiologischen 
elektrophysiologischen Prozesses. Richtungen der momentanen elektromagnetischen Kraft (EMF) und der 
Stromdichte, elektrische Atria und ventrikulare Aktivitat. 



Schlusse hinsichtlich verschiedener Fomien von myocardialer Ischemie. supraventrikularer und ventrikularer 
RhythmusstOrungsrisiken, Rate des Herzversagens, Veranderungen des elektrophysiologischen Prozesses 
unter Behandlung Oder diagnostische Tests. 



Tabelle 1 
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Patentanspruche 

1 . Computerbasiertes Verfahren zur automatischen Aufbereitung von Daten bio- 
magnetischer Felder, insbesondere von nnagnetokardiographischen Daten, umfassend 

5 die folgenden, auf einer OberflSchendichte magnetischer Momente (Schicht magnetic 
scher Dipole) Oder, was Equivalent hierzu ist, einer Funktion der Strdme als physikali- 
sches und mathematisches Modell der Quellen eines biomagnetischen Feldes beru- 
henden Schritte: 

Aufstellen einer die Oberflachendichte der magnetischen Momente betreffen- 
10 den Integralgleichung, deren rechtes Glied die von einem Gradiometer gemessene 
zweite Ableitung der magnetischen Feldinduktion in Normalenrichtung zur Mefiebene 
(a^Bz/az^) darstellt; 

Bestimmen analytischer Ausdrucke fur Faktoren einer Matrix A, die den integra- 
len Operator dier genannten Intergralgleichung anndhert, und Bereclinung dieser Ma- 
15 trix; 

Interpolation der gemessenen Werte der Funktion y^S^Bz/Sz^ in den Knoten 
eines vorzugsweise klein dimensionierten Gitters; 

Losung gemaB Tickhonov eines iinearen algebraischen Gleichungssystems 
Ax = y, wobei x die Oberflachendichte der magnetischen Momente ist; 
20 - Aufbau einer Niveaulinienkarte der Oberflachendichte der magnetischen Mo- 
mente Oder, was Equivalent hierzu ist, einer Stromlinienkarte und Auslesen der Karte in 
eine Speichereinheit oder eine Ausgabeeinheit. 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, da& der Gradienten der 
25 Oberflachendichte der magnetischen Momente automatisch berechnet und eine Karte 

der Stromdichten aufgebaut und in eine Speichereinheit Oder eine Ausgabeeinheit 
ausgelesen wird. 

3. Computerbasiertes Verfahren zur automatischen Aufbereitung von Daten bio- 
30 magnetischer Felder, insbesondere von magnetokardiographischen Daten, insbeson- 
dere nach Anspruch 1, dadurch gekennzeichnet, da& zumindest ausgewShlte Daten 
nach vorgegebenen Kriterien automatisch mit wenigstens einem vordefinierten Nor- 
malwert verglichen werden und daB bei Abweichung der Daten von dem Normalwert 
um ein vorgegebenes Ma& ein auf einer Ausgabeeinrichtung ausgebbares die Abwei- 

35 Chung signalisierendes Signal erzeugt wird. 
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4. Verfahren nach Anspruch 3. dadurch gekennzeichnet, dali Gr6(ie der Felder 
unter der P-Welle und dem QRS-Komplex berechnet werden. Die GroBe dieser Felder 
spiegelt die Energie wider, die wahrend der Anregung der Vorkammem und Ventriculi 
des Herzens erzeugt wird. Femer kann das Verh§ltnis dieser Felder berechnet werden. 

5 urn so die elektrische Aktlvitat der Vorkammem im Vergleich zu elektrischen Aktlvitat 
der Ventrikel abzuschatzen. Es hat sich gezeigt. da(S dieser Parameter mit dem Grad 
des Herzversagens in Beziehung steht. 

5. Verfahren nach einem der Anspriiche 1 bis 4. dadurch gekennzeichnet. da(i 
10 einer der in der Beschreibung erwahnten Parameter aus den erfaliten Daten automa- 

tisch berechnet und mit einem Normalwert verglichen wird. 

6. Verfahren nach einem der Anspriiche 1 bis 5, dadurch gekennzeichnet, da(i die 
erfaliten Daten automatisch unterschiedliche Analysestufen. wobei die Komplexitat der 

15 Analyse auf jeder Stufe zunimmt. durchlaufen. 

7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, daU wenigstens eine der 
folgenden Analysestufen durchlaufen wird: Amplituden-Zeit-Analyse verschiedener 
MCG-Kurven; Analyse hoch-aufl6sender MCG; Analyse summierter MCG und Vek- 

20 tor-MCG (VMCG); Qualitative Analyse der magnetischen Feld-Verteilungskarten; 

Quantitative Analyse der magnetischen Feld-Verteilungskarteri; Analyse der effektiven 
Dipol-Charakteristiken; Analyse der zweidimensionalen Ladungsverteilung; Analyse 
der dreidimensionalen Ladungsverteilung. 

25 8. Computerbasiertes Verfahren zur automatischen Aufbereitung magnetokardio- 
graphischer Daten, umfassend die Schritte: 

Abschdtzung der wechselseitigen Anordnung der maximalen positiven und ne- 
gativen Extreme in jeder magnetischen Feldverteilungskarte. 

Abschdtzen der Richtung des dquivalenten Stromdipols in jeder Karte unter 
30 Verwendung der Rechte-Hand-Regel. 

Abschatzen der Anwesenheit zusStzlicher Extrema in jeder magnetischen Feld- 
verteilungskarte und Dauer deren Existenz wShrend des Repolarisationsprozesses. 

Berechnung des Verhaitnisses der maximalen negativen zu den maximalen 
positiven Extrema in jeder Karte und der Standardableitung dieser Parameter wahrend 
35 des Verfahrens. 
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Abschatzung der effektiven Dipoltiefenkurve wahrend der Repolarisation, 
Abschatzung der Verteilung und wechselseitigen Anordnung einer Stromlinie 
und der Richtungen der Stromdichtevectoren in jeder Karte wahrend der Repolarisati- 
on. 

5 - Abschatzung der Kurven durchschnittlicher und maximaler Stronndichte wah- 
rend der Repolarisation. 

automatisches Bewerten der ventrikularen Repolarisation hinsichtlich der Kon- 
fonnitat des aktuellen Repolarisationsprozesses im Vergleich zu nornnalen Parametem. 

10 9. Computerbasiertes Verfahren zur automatischen Aufbereitung nnagnetokardio- 
graphischer Daten, umfassend die folgenden Schrltte 

die wechelseitige Anordnung der nfiaximalen positiven und negativen Extreme 
In jeder magnetischen Feldverteilungskarte wird abgeschatzt. 

GemaU der Rechte-Hand-Regel wird die Richtung des aquivalenten Stromdi- 
15 pols in jeder Karte abgeschatzt. 

Das Vorhandensein zusatzlicher Extrema in jeder magnetischen Feldvertei- 
lungskarte und die Dauer deren Vorhandensein wird wahrend des Depolarisationspro- 
zesses abgeschatzt. 

Die effektive Dipoltiefenkurve wahrend der Depolarisation wird abgeschatzt. 
20 - Die Verteilung und wechselseitige Anordnung von Stromlinien und die Richtun- 
gen der Stromdichtevectoren in jeder Karte wahrend der Depolarisation wird abge- 
schatzt. 

Die Dauer der Existenz zusatzlicher Strombezirke wahrend des Depolarisati- 
onsprozesses wird abgeschatzt. 
25 - Die Kurven der durchschnittlichen und maximalen Stromdichte wahrend der 
Depolarisation wird abgeschatzt. 

Automatisches Bewerten der Konformitat des aktuellen Depolarisationsprozes- 
ses mit nomialen Parametem. 

30 10. Computerbasiertes Verfahren zur automatischen Aufbereitung magnetokardio- 
graphischer Daten, umfassend die folgenden Schritte 

die wechselseitige Anordnung der maximalen positiven und negativen Extrema 
. in jeder magnetischen Feldverteilungskarte wird abgeschatzt. 

GemaS der Rechte-Hand-Regel wird die Richtung der aquivalenten Stromdipole 
35 in jeder Karte abgeschatzt. 
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Die Anwesenheit zus^tzlicher Extrema in jeder magnetischen Feldverteilungs- 
karte und die Dauer deren Existenz wird wahrend des Repolarisationsprozesses abge- 
schatzt. 

Die Ordnungszahl der Karte beginnend vom Zeitpunkt. zu dem die Richtung 
5 des aquivalenten Stromdipols der ventrikularen Repolarisation stabil normal wurde, 
d.h. nach links unten gerichtet. 

Die effektive Dipoltiefe im J-Punkt wird abgeschStzt. 

Die Verteilung und wechselseitige Anordnung der Stromlinien und der Richtun- 
gen der Stromdichtevectoren in jeder Karte wShrend der Repolarisation wird abge- 
10 schatzt. 

Die Ordnungszahl der Karte beginnend vom Zeitpunkt, zu dem die Richtung der 
maximalen Stromdichtevectoren der ventrikularen Repolarisation stabil normal wurde, 
d,h. nach links unten gerichtet, wird bestimmt. 

Die Dauer der Existenz zusStzlicher Strombereiche wShrend des Repolarisatl- 
1 5 onsprozesses wird abgeschatzt. 

Die Kurven der durchschnittlichen und der maximalen Stromdichte wahrend der 
Repolarisation am Beginn des ST-T-lntervalls werden abgeschatzt 

Die Verhaitnisse der maximalen und der durchschnittlichen Stromdichten im 
QRS-Maximum zu dem im T-Maximum werden abgeschatzt. 
20 Automatisches Bewerten der Abschatrzungen hinsichtlich des Vorhandenseins 

und der Schwere einer myocardialen Ischemie, insbesondere bei bei Patienten mit 
nicht-verandertem Ruhe-EKG umfassend die foigenden Schritte: 

1 1 . Computerbasiertes Verfahren zur automatischen Aufbereitung magnetokardio- 
25 graphischer Oaten, umfassend die foigenden Schritte 

die wechselseitige Anordnung der maximalen positiven und negativen Extrema 
In jeder magnetischen Feldverteilungskarte wird abgeschatzt. 

Gemaii der Rechte-Hand-Regel wird die Richtung der aquivalenten Stromdipole 
in jeder Karte abgeschatzt. 
30 - Emnitteln der Zeitperioden, wahrend deren die Richtung der effektiven Dipole 1 . 
nach rechts Oder nach rechts unten, 2. nach links unten. 3. nach rechts Oder Oder nach 
oben gerichtet blelbt. 

Die Anwesenheit zusatzlicher Extrema in jeder magnetischen Feldverteilungs- 
karte und die Dauer deren Existenz wahrend des Depolarisationsprozesses wird abge- 
35 schatzt. 
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Die relative Amplitude der nach oben und nach unten gerichteten Bewegung 
der effektiven Oipole wdhrend der Depolarisation wird abgeschdtzt. 

Die Verteilung und wechselseitige Anordnung einer Stromlinie und die Richtun- 
gen der Stromdichtevectoren in jeder Karle wShrend der Depolarisation wird abge- 
5 schdtzt. 

Ermitteln der Zeitperiode, wShrend welcher die Richtung der maximalen 
Stromdichtevectoren 1 . nacfi rechts Oder nach rechts unten gerichtet, 2. nach links 
unten gerichtet, drittens nach rechts Oder nach rechts oben gerichtet blelbt. 

Abschatzung der Dauer der Existenz zusStzlicher Strombereiche wahrend des 
10 Depolarisationsprozesses wird abgeschatzt. 

Die Spitzenwerte der durchschnittlichen und maximalen Stromdichte wahrend 
der Depolarisation und die Form der entsprechenden Kurven wird abgeschatzt. 

Automatisches Bewerten der Abschatzungen hinsichtlich des Vorhandenseins 
myocardialer Nekrose. ^ 

15 

12. Computerbasiertes Verfahren zur automatischen Aufbereitung von Daten bio- 
magnetischer Felder, insbesondere von magnetokardiographischen Daten, wie aus der 
Beschreibung ersichtlich. 

20 13. Verfahren nach einem der AnsprQche 1 bis 12, dadurch gekennzeichnet, da& 
belm Vergleich der betrachteten Daten mit den Normalwerten die betrachteten Daten 
automatlsch in einer selbstlernenden Datenbank gespeichert werden und nach vor- 
gebbaren Kriterien zur fortlaufenden Bildung und Oberpriifung von Normalwerten und 
tolerablen Abweichungen der Normalwerte hinzugezogen werden kdnnen. 

25 

14. Anwendung eines der beanspruchten und/oder beschriebenen Verfahren zur 
automatischen Analyse der Daten eines Patentienten, vorzugsweise der Daten grO&ere 
Patentientenmengen im Rahmen, insbesondere eines Screeningtest, zur Evaluation 
des Risikos, an bestimmten Krankheiten, insbesondere einer der in der Anmeldung 
30 genannten Herzkrankhelten zu erkranken. 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



International application No. 

PCT/DEOO/02930 



I. Basis of the report 



1. With regard to the elements of the international application:* 
I I the international application as originally filed 

1^ the description: 

pages [ 1-22 

pages , 

pages 



, as originally filed 

, filed with the demand 



, filed with the letter of 



the claims: 

pages 

pages 

pages 

pages 



M4 



, as originally filed 

, as amended (together with any statement under Article 19 
, filed with the demand 



filed with the letter of 



the drawings: 

pages 

pages 

pages 



, as originally filed 

, filed with the demand 



filed with the letter of 



I I the sequence listing part of the description: 

pages 

pages 

pages 



, as originally filed 

, filed with the demand 



filed with the letter of 



2. With regard to the language, all the elements marked above were available or furnished to this Authority in the language in which 
the international application was filed, unless otherwise indicated under this item. 

These elements were available or furnished to this Authority in the following language which is: 

O the language of a translation furnished for the purposes of international search (under Rule 23. 1 (b)). 
n the language of publication of the international application (under Rule 48.3(b)). 

Q the language of the translation furnished for the purposes of international preliminary examination (under Rule 55.2 and/ 
or 55.3). 

3. With regard to any nucleotide and/or amino acid sequence disclosed in the international application, the international 
preliminary examination was carried out on the basis of the sequence listing: 

n contained in the international application in written form. 

Q filed together with the international application in computer readable form. 

furnished subsequently to this Authority in written form. 

0 furnished subsequently to this Authority in computer readable form. 

n The statement that the subsequently furnished written sequence listing does not go beyond the disclosure in the 

international application as filed has been furnished. 

n The statement that the information recorded in computer readable form is identical to the written sequence listing has 
been furnished. 

n The amendments have resulted in the cancellation of: 

I I the description, pages , 

I I the claims, Nos. 

I I the drawings, sheets/fig , 

5 rn I'his report has been established as if (some of) the amendments had not been made, since they have been considered to go 
— ' beyond the disclosure as filed, as indicated in the Supplemental Box (Rule 70.2(c)).** 

♦ Replacement sheets which have been furnished to the receiving Office in response to an invitation under Article 14 are referred to 
in this report as "originally filed" and are not annexed to this report since they do not contain ametuiments (Rule 70.16 
and 70 J 7). 

^*Any replacement sheet containing such amendments must be referred to under item J and annexed to this report. 
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III. Non-establishment of opinion with regard to novelty, inventive step and industrial applicability 



1. The questions wliether the claimed invention appears to be novel, to involve an inventive step (to be non obvious) or to be 
industrially applicable have not been examined in respect of: 

I I the entire international application, 

^ claims Nos. N14 



because: 



□ the said international application, or the said claims Nos. 
relate to the following subject matter which does not requT 



require an international preliminary examination (specify): 



^ the description, claims or drawings (indicate particular elements below) or said claims Nos. M4 

^ are so unclear that no meanitigful opinion could be formed (specify): 

See supplemental Box 



□ the claims, or said claims Nos. ^ are so inadequately supported 
by the description that no meaningful opinion could be formed. 

I I no international search report has been established for said claims Nos. 



2. A meaningful international preliminary examination cannot be carried out due to the failure of the nucleotide and/or amino acid 
sequence listing to comply with the standard provided for in Annex C of the Administrative Instructions: 



I ] the written form has not been furnished or does not comply with the standard. 

[~[ the computer readable form has not been furnished or does not comply with the standard. 
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Supplemental Box 

(To be used when the space in any of the preceding boxes is not sufficient) 

Continuation of: Box III.l 

The following document is cited by the applicant in the 
description (the reference niamber "Dl" will be retained 
throughout the remainder of the proceedings) : 

Dl : Medical & Biological Engineering & Computing, Vol. 33, 
No. 5, September 1995, GB, pages 543-651, R. Killmann 
et al.: "Localisation of myocardial ischaemia from the 
magnetocardiogram using current density reconstruction 
method: computer simulation study" 

1. In view of the objections raised in Box VIII. below, the 

international preliminary examining authority is unable to 
establish an expert opinion under the Regulations (PCT 
Article 34(4) (a) (ii) ) with regard to novelty, inventive 
step and industrial applicability for the subject matter 
of Claims 1-14. 

1.1 However, if Claims 1-14 were amended as suggested in Box 
VIII, they would, insofar as can be understood, meet the 
requirements of novelty, inventive step and industrial 
applicability according to PCT Article 33(2), (3) and (4). 
The reasons for this are as follows: 

1.2 The invention relates to a computer-based method for the 
automatic processing of biomagnetic field data, in 
particular magnetocardiographic data. 

Specifically the invention relates to a method which uses 
one or more detectors, in particular a SQUID detector, to 
allow biomagnetic fields measured in a plane or at several 
points in a plane to be automatically visualised, for 
example in the form of magnetic field maps, and to allow 
the information contained in the field to be automatically 
analysed and classified. 
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Supplemental Box 

(To be used when the space in any of the preceding boxes is not sufficient) 
Continuation of: Box 1 1 1 . 1 



The development of highly sensitive magnetographs , 
particularly those which can be used in unshielded 
environments, requires new methods for processing the 
measured data, since the methods known to date are 
extremely time-consuming and there is only a small number 
of specialists able to carry them out. As far as possible, 
the data should be processed automatically so that the 
diagnostically relevant information in the measured 
magnetic fields can be easily identified and interpreted by 
doctors. In a preferred embodiment, the method should 
compare the diagnostically relevant information contained 
in the data with stored information and automatically 
classify it in order to assist the doctor in his diagnosis. 

Certain information needs to be visualised; for example, 
in the form of current line maps. 

If it is desired to show biomagnetic currents, especially 
those in the heart, which have proven particularly relevant 
in clinical diagnostics, an "inverse problem" needs to be 
solved; that is, the magnetic fields generated by the 
biomagnetic currents must be measured and used to calculate 
the currents which generate the fields . 

Given the complexity of the inverse problem in the 
localisation of magnetic field sources, the solution makes 
various simplifying assumptions. One of these assumptions, 
for example, is that all the sources are positioned in the 
same plane and at a particular distance from a measuring 
surface. The results obtained depend on which physical 
model and hence which mathematical model is chosen as the 
model of the magnetic field sources. If a flat source 
region is chosen, the model chosen for the sources is 
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Supplemental Box 

(To be used when the space in any of the preceding boxes is not sufficient) 

Continuation of: Box III.l 

normally the distribution of current densities on a 
surface, which represents the measured vector quantity. 
Solving the inverse problem produces a map of directional 
arrows, the magnitudes of which are proportional to the 
amount of a current density. This type of current imaging 
causes certain difficulties. On the side of the calculated 
region where the modulo value of the current densities is 
several times smaller than the maximum it is not readily 
possible to estimate the magnitudes and directions of the 
arrows on the arrow map. The low-intensity sources are 
> therefore not taken into account in the analysis. 

A method for automatically visualising cardiac currents 
using current lines in a plane is described in document Dl. 
In the method according to Dl, the solution to the inverse 
problem is a map of directional arrows as indicated above, 
in which the magnitudes of the arrows are proportional to 
the modulus of the current density. 

The disadvantage of the cardiac current imaging method 
known from Dl using current lines in a plane is that in 
regions in which the current densities are modulus 3 or 
more below the maximum it is not readily possible to make a 
visual assessment to determine magnitude and direction in 
the arrow pattern generated by the method. Experience has 
shown that to ensure a reasonable error estimate it is only 
possible to visualise one or two particularly intense 
regions of a magnetic field. The other, less intensive 
regions cannot be perceived visually. 

The problem addressed by the invention is therefore that of 
providing a method for processing the data from a measured 
biomagnetic field, whereby the currents generating the 
biomagnetic field can be visualised by way of current lines 
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Supplemental Box 

(T o be used when the space in any of the preceding boxes is not suflQcient) 

Continuation of: Box III.l 

in the plane so that visual analyses can be readily made by 
doctors or trained persons . 

The problem is solved by a method for automatically 
processing biomagnetic field data, in particular 
magnetocardiographic data, comprising the following steps 
involving a surface density of magnetic moments: 

(i) establishing an integral equation for the surface 
density of the magnetic moments, in which the 
right-hand side represents the second derivative 
of the magnetic field induction in the direction 
normal to the measuring plane as measured by a 
gradiometer (d^Bjdz^) ; 

(ii) finding analytic expressions for factors of a 
matrix A which approximates the integral operator of 
said integral equation, and calculating said matrix; 

(iii) interpolating the measured values for the function 
y = a^Bj/az^ in the node of a preferably small grid; 

(iv) using the Tikhonov method to solve a system of linear 
algebraic equations Ax = y, where x is the surface 
density of the magnetic moments; 

(v) constructing a level line map of the surface density 
of the magnetic moments. 

As far as can be understood, the solution proposed by the 
present application for the aforementioned problem appears 
to be novel, and for the reasons given above it also appears 
to involve an inventive step. 
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VII. Certain defects in the international application 

The following defects in the form or contents of the international application have been noted: 

2. The claims include references to the description and 

drawings. Under PCT Rule 6.2(a), claims may only include 
references when they are absolutely necessary^ which is 
not the case here. 

2.1 The features of the claims are not followed by reference 
signs in parentheses (PCT Rule 6.2(b)). 
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VIII. Certain observations on the international application 

The following observations on the clarity of the claims, description, and drawings or on the question whether the claims are fully 
supported by the description, are made: 

3. Claims 1, 8, 9, 10, 11, 12 and 14 have been drafted as 

separate independent claims, yet they actually appear to 
relate to the same subject matter, the only apparent 
difference being in the definition of the subject matter 
for which protection is sought and in the terminology used 
to define the features of the claimed subject matter. The 
claims therefore lack concision. Moreover, the claims 
display an overall lack of clarity because the large 
number of independent claims makes it hard, if not 
impossible, to identify the subject matter for which 
protection is sought, and it is therefore unreasonably 
difficult for a third party to determine the scope of 
protection. Claims 1, 8, 9, 10, 11, 12 and 14 therefore 
fail to meet the requirements of PCT Article 6. 

3.1 This objection could be overcome by submitting an 

amended set of claims in which the subject matter is 
defined by a single independent claim, supplemented if 
necessary by dependent claims covering the optional 
features (PCT Rule 6.4). 
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WEITERES VORGEHEN 



siehe Mitteilung uber die Ubersendung des Intemationalen 
vortSufigen Prufungsberichts (Formblatt PCT/I PEA/41 6) 



Internationales Aktenzelchen 
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Internationales Anmeldedatum(T ag/Monat/Jahr) 
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Prioritatsdatum (Tag/Monat/Tag) 
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Intematlonate Patentklassifikation (IPK) Oder nationale Klasslfikation und IPK 
G06F17/00 



Anmelder 

SQUID AG et al. 



1. DIeser internationale vorlaufige Prufungsbericht wurde von der mit der intemationalen vorlSuflgen Priifung beauftragten 
Behorde erstellt und wird denn Anmelder gemaB Artikel 36 Qbermittelt. 

2. Dieser BERICHT umfaBt insgesamt 9 Blotter eInschlieBlich dieses Deckblatts. 

□ AuBerdem liegen dem Berlcht ANLAGEN bel; dabei handelt es sich urn Blatter mit Beschrelbungen, Anspruchen 
und/oder Zeichnungen, die geandert wurden und diesem Bericht zugrunde liegen, und/oder Blatter mit vor dieser 
Behorde vorgenommenen Berichtigungen (siehe Regel 70.16 und Abschnitt 607 der Verwaltungsrichtlinien zum PCT). 

Diese Aniagen umfassen Insgesamt Blatter. 



3. Dieser Bericht enthalt Angaben zu folgenden Punkten: 

I ^ Grundlage des Berichts 
Priorltat 

Kelne Erstellung eines Gutachtens uber Neuhelt. erfinderische Tatigkeit und gewerbliche Anwendbarkelt 
Mangelnde EInheitlichkelt der Erfindung 

Begrundete Feststellung nach Artikel 35(2) hinsichtlich der Neuheit, der erfinderischen Tatigkeit und der 
gewerblichen Anwendbarkeit; Unterlagen und Erklarungen zur Stutzung dieser Feststellung 



II 
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IV 
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Formblatt PCT/IPEA/409 (Deckblatt) (Januar 1994) 




INTERNATIONALER VORLAUFIGER 

PRUFUNGSBERICHT 



Internationales Aktenzeichen PCT/DEOO/02930 



. I. Grundlage des Berichts 

1 . Hinsichtllch der Bestandteile der internationalen Anmeldung {Ersatzblatter, die dem Anmeldeamt auf eine 
Aufforderung nach Artikel 14 bin vorgelegt warden, gelten im Rahmen dieses Benefits als "ursprunglich 
eingereicht" und sind ihm nicht beigefugt, weil sie keine Anderungen enthalten (Regein 70. 16 and 70. 17)): 
Beschreibung, Seiten: 

1-22 ursprungliche Fassung 

Patentanspruche, Nr.: 

1 -1 4 ursprungliche Fassung 

Zeichnungen, Blatter: 

1/1 ursprungliche Fassung 



2. Hinsichtllch der Sprache: Alle vorstehend genannten Bestandteile standen der Behdrde in der Sprache, in der 
die internationale Anmeldung eingereicht worden ist, zur Verfugung Oder wurden in dieser eingereicht, sofern 
unter diesem Punkt nichts anderes angegeben ist. 

Die Bestandteile standen der Behorde in der Sprache: zur Verfugung bzw. wurden in dieser Sprache 
eingereicht; dabei handelt es sich urn 

□ die Sprache der Ubersetzung, die fur die Zwecke der internationalen Recherche eingereicht worden ist (nach 

Regel 23.1(b)). 

□ die Veroffentlichungssprache der internationalen Annneldung (nach Regel 48.3(b)). 

□ die Sprache der Ubersetzung, die fur die Zwecke der internationalen vorlaufigen Prufung eingereicht worden 
ist (nach Regel 55.2 und/oder 55.3). 

3. Hinsichtllch der in der internationalen Anmeldung offenbarten Nucleotid- und/oder Aminosauresequenz ist die 
internationale vorlaufige Prufung auf der Grundlage des Sequenzprotokoils durchgefuhrt worden, das: 

□ in der internationalen Anmeldung In schriftlicher Form enthalten ist. 

□ zusammen mit der internationalen Anmeldung in computerlesbarer Form eingereicht worden ist. 

□ bei der Behdrde nachtragllch in schriftlicher Form eingereicht worden ist. 

□ bei der Behorde nachtragllch in computerlesbarer Form eingereicht worden ist. 

□ Die Erklarung, daB das nachtragllch eingereichte schriftliche Sequenzprotokoll nicht uber den 
Offenbarungsgehalt der internationalen Anmeldung im Anmeldezeitpunkt hinausgeht, wurde vorgelegt. 

□ Die Erklarung, daB die in computerlesbarer Form erfassten Informationen dem schriftlichen 
Sequenzprotokoll entsprechen, wurde vorgelegt. 

4. Aufgrund der Anderungen sind folgende Unterlagen fortgefallen: 
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□ 



Beschreibung, 

Anspruche, 

Zeichnungen, 



Seiten: 



□ 



Nr.: 



□ 



Blatt: 



5. □ Dieser Bericht 1st ohne Berucksichtigung (von einigen) der Anderungen erstellt worden, da diese aus den 

angegebenen Grunden nach Auffassung der Behorde uber den Offenbarungsgehalt In der ursprunglich 
eingereichten Fassung hinausgehen (Regel 70.2(c)). 

(Auf Ersatzblatter, die solche Anderungen enthaiten, ist unter Punkt 1 hinzuweisen;sie sind diesem Bericht 
belzufugen), 

6. Etwaige zusatzliche Bemerkungen: 

III. Kelne Erstellung eines Gutachtens iiber Neuheit, erfinderische Tatigkeit und gewerbliche Anwendbarkeit 

1 . Folgende Telle der Anmeldung wurden nicht daraufhin gepruft, ob die beanspruchte Erfindung als neu, auf 
erfinderischer Tatigkeit beruhend (nicht offensichtlich) und gewerbllch anwendbar anzusehen ist: 

□ die gesamte Internationale Anmeldung. 
la Anspruche Nr. 1-14. 



□ Die gesamte internationale Anmeldung, bzw. die obengenannten Anspruche Nr. bezlehen sich auf den 
nachstehenden Gegenstand. fur den kelne internationale vorlaufige Prufung durchgefuhrt werden braucht 
(genaue Angaben): 



la Die Beschreibung, die Anspruche oder die Zeichnungen (machen Sle hierzu nachstehend genaue Angaben) 
Oder die obengenannten Anspruche Nr. 1-14 sind so unklar, daB kein sinnvolles Gutachten erstellt werden 
konnte (genaue Angaben): 
siehe Belblatt 

□ Die Anspruche bzw. die obengenannten Anspruche Nr. sind so unzureichend durch die Beschreibung 
gestutzt, daB kein sinnvolles Gutachten erstellt werden konnte. 

□ Fur die obengenannten Anspruche Nr. wurde kein internationaler Recherchenbericht erstellt. 

2. Eine sinnvolle internationale vorlaufige Prufung kann nicht durchgefuhrt werden, weil das Protokoll der Nukleotld- 
und/oder Aminosauresequenzen nicht dem in Aniage C der Verwaltungsvorschriften vorgeschriebenen Standard 



Begrundung: 



\ 



entspricht: 



□ Die schriftliche Form wurde nicht eingereicht bzw. entspricht nicht dem Standard. 

□ Die computerlesbare Form wurde nicht eingereicht bzw. entspricht nicht dem Standard. 
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Vli. Bestimmte Mangel der mternatlonalen Anmeldung 

Es wurde festgestellt. da3 die intematlonale Anmeldung nach Fomn Oder Inhalt folgende Mangel aufweist: 
siehe Beiblatt 



VIII. Bestimmte Bemerkungen zur internationalen Anmeldung 

Zur Klarheit der Patentanspruche, der Beschreibung und der Zeichnungen oder zu der Frage, ob die Anspriiche 
in vollem Umfang durch die Besclireibung gestutzt werden, 1st folgendes zu bemerken: 
siehe Beibiatt 
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Das folgende Dokument (D) wurde von der Anmelderin in der Beschreibung 
zitlert. Die Numerierung wird auch im weiteren Verfahren beibehalten. 

D1 : Medical & Biological Engineering & Computing, Vol. 33 N°5 September 
1995, GB, pages 643-651, R.Killmann et al: "Localisation of myocardial 
Ischaemia from the magnetocardlogram using current density reconstruction 
method: computer simulation study" 

Zu Punkt III 

Keine Erstellung eines Gutachtens uber Neuhelt, erfinderlsche Tatlgkeit und 
gewerbliche Anwendbarkelt 

Wegen der Punkt 8 unter erhobenen Einwande, ist es fiir die mit der 
intemationalen vorlaufigen Prufung beauftragte Behorde nicht mogllch, nach der 
Ausfuhrungsordnung (siehe Artikel 34 (4) (a) (ii) PCT), ein Gutachten iiber 
Neuheit, erflnderische Tatigkeit und gewerbliche Anwendbarkelt in bezug auf den 
Gegenstand der Anspriiche 1 bis 14 zu erstellen. 

Sollten Anspruche 1 bis 14 jedoch wie in Punkt 8 vorgeschlagen geandert 
werden, wiirden sie, soweit zu verstehen ist, die Erfordernisse der Neuheit, der 
erfinderischen Tatigkeit und gewerbliche Anwendbarkeit nach Artikel 33 (2), (3) 
und (4) PCT erfullen; die Grunde dafur sind die folgenden. 

1 .2 Die Erfindung betrifft ein computerbasiertes Verfahren zur automatischen 
Aufbereitung von Daten biomagnetischer Felder, insbesondere von 
magnetokardiographischen Daten. 

Insbesondere betrifft die Erfindung ein Verfahren, das es eriaubt, daf3 mittels 
eines Oder mehrerer Detektoren, insbesondere mittels eines SQUID-Detektors in 
einer Ebene bzw. an mehreren Punkten einer Ebene aufgenommene 
biomagnetische Felder automatisch z. B. in Form von Magnetfeldkarten zu 
visualisieren und die in dem Feld enthaltenen Informationen automatisch zu 
evaluieren und zu klassifizieren. 

Die Entwicklung hoch-sensitiver Magnetographen, insbesondere solcher, die in 
nicht-abgeschirmten Umgebungen eingesetzt werden konnen, erfordert neue 



1 



1.1 
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Verfahren zur Aufbereitung der erfaBten Daten, da die bislang bekannten 
Verfahren auBerst zeitaufwendig sind und nur von wenigen Spezialisten 
ausgefiihrt werden konnen. Die Daten sollen weitestgeliend automatisch so 
aufbereitet werden, daB die in den erfaBten Magnetfeldem enthaltenen 
diagnostisch relevanten Informationen von einem Arzt leicht erkannt und bewertet 
werden konnen. Bei einer bevorzugten Ausfiihrungsform soli das Verfahren die in 
den Daten enthaltenen diagnostisch relevanten Informationen mit gespeicherten 
Informationen vergleichen und automatisch klassifizieren, um so den Arzt bei 
seiner Diagnose zu unterstutzen. 

Bestimmte Informationen sollen z. B. in Form von Stromlinienkarten visualisiert 
werden. 

Sollen biomagnetische Strome, insbesondere die Strome im Herzen, die in der 
klinischen Diagnostik erwiesenermaBen eine besonders hohe Relevanz besitzen, 
dargestellt werden, muB ein sogenanntes "inverses Problem" gelost werden, das 
heiBt, daB aus den erfaBten Magnetfeldem, die von den biomagnetischen 
Stromen hervorgerufen werden, die die Felder hervorrufenden Strome berechnet 
werden miissen. 

In Anbetracht der Komplexitat eines inversen Problems bei der Lokalisation 
magnetischer Feldquellen werden verschiedene vereinfachende Annahmen bei 
der Losung gemacht, Eine dieser Annahmen ist beispielsweise, daB alle Quellen 
in derselben Ebene positioniert und in einer bestimmten Distanz von einer 
MeBoberflache angeordnet sind. Die erzielten Resultate davon ab, welches 
physikalische Modell und entsprechend welches mathematische Modell man als 
Modell der magnetischen Feld- quellen wahlt. Falls man einen flachen 
Quellenbereich wahlt, wird ublichenA^eise die Verteilung der Stromdichten auf 
einer Oberflache als Modell der Quellen gewahit, die die erfaBte VektorgroBe 
darstellt. Nach Losung eines inversen Problems erhalt man eine Karte von 
gerichteten Pfeilen, deren GroBe proportional zum Betrag einer Stromdichte ist. 
Das Abbilden von Stromen in einer solchen Weise bedingt bestimmte 
Schwierigkeiten. Auf Seiten des berechneten Bereiches, wo der Modulo-Wert der 
Stromdichten mehrere Male kleiner ist als das Maximum, ist es nicht ohne 
weiteres moglich, auf der Karte von Pfeilen die GroBe und Richtung der Pfeile 
abzuschatzen. Die Quellen mit kleiner Intensitat fallen so aus der Untersuchung. 
Ein Verfahren zur automatischen Visualisieren von Herzstromen durch 
Stromlinien in der Ebene wird in D1 beschrieben. GemaB dem in D1 
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beschriebenen Verfahren erhalt man als Resultat der Losung des inversen 
Problems eine der genannten Karten von orientierten Pfeilen, deren GroBe 
proportional zum Modulus der Stromdichte ist. 

Das aus dem D1 bekannte Verfahren zum Abbilden von Herzstromen durch 
Stromlinien in der Ebene hat jedoch den Nachteil. daB in Bereichen, in denen die 
Stromdichten Modulus 3 und mehr geringer sind als das Maximum, eine visuelle 
Beurteilung des mil dem Verfahren erzeugten Musters von Pfeilen hinsichtlich 
ihrer GroBe und Richtung nicht mehr ohne weiteres moglich ist. In der Praxis hat 
sich gezeigt, daB im Hinblick auf eine verniinftige Fehlerabschatzung lediglich ein 
bis zwei besonders intensive Bereiche eine magnetischen Feldes visualisiert 
werden konnen. 

Die ubrigen, weniger intensiven Bereiche konnen nicht sichtbar wahrgenommen 
werden. 

Davon ausgehend liegt der Erfindung zum einen die Aufgabe zugrunde, ein 
Verfahren zum Aufbereiten der Daten des erfaBten biomagnetischen Feldes 
anzugeben, mittels welchem die das biomagnetische Feid hervorrufenden Strome 
durch Stromlinien in der Ebene visualisiert werden konnen und so einer visuellen 
Auswertung durch geschultes Personal Oder einen Arzt in besonders einfacher 
Weise zuganglich werden. 

Die Aufgabe wird gelost von einem Verfahren zur automatischen Aufbereitung von 
Daten biomagnetischer Felder, insbesondere von magnetokardiographischen 
Daten, umfassend die folgenden, auf einer Oberflachendichte magnetischer 
Momente, beruhenden Schritte: 

i) Aufstellen einer die Oberflachendichte der magnetischen Momente 
betreffenden Integralgleichung, deren rechtes Glied die von einem 
Gradiometer gemessene zweite Ableitung der magnetischen Feldinduktion 
in Normalenrichtung zur MeBebene, d^BJdz^, darstellt ; 

ii) Bestimmen analytischer Ausdriicke fur Faktoren einer Matrix A, die den 
integralen Operator der genannten Integralgleichung annahert, und 
Berechnung dieser Matrix ; 

iii) Interpolation der gemessenen Werte der Funktion y=d^BJdz^ in den Knoten 
eines vorzugsweise klein dimensionierten Gitters; 

iv) Losung gemaB Tickhonov eines linearen algebraischen Gleichungssystems 
Ax = y, wobei x die Oberflachendichte der magnetischen Momente ist; 
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v) Aufbau einer Niveaulinienkarte der Oberflachendichte der magnetlschen 
Momente. 

Soweit zu verstehen ist, scheint es, daB die fiir die obengenannte Aufgabe 
vorgeschlagene Losung der vorliegenden Anmeldung somit neu ist und, aus den 
vorherigen Grunden, auf einer erfinderisclien Tatigl<eit beruht. 

Zu Punkt VII 

Bestimmte Mangel der internationalen Anmeldung 

Die Anspruche entlialten Bezugnahmen auf die Beschreibung bzw. die Zeichnun- 
gen. GemaB Regal 6.2 a) PCT diirfen Anspruche nur dann Bezugnahmen 
enthalten, wenn dies unbedingt erforderlich ist, was hier nicht der Fail ist. 

Die Merkmale der Anspruche sind nicht mit in Klammern gesetzten 
Bezugszeichen versehen worden (Regel 6.2 b) PCT). 



Zu Punkt Vlil 

Bestimmte Bemerkungen zur internationalen Anmeldung 

Die Anspruche 1, 8, 9, 10, 1 1, 12 und 14 wurden zwar als getrennte, unabhangige 
Anspruche abgefaBt, sie scheinen sich aber tatsachlich auf ein und denselben 
Gegenstand zu beziehen und unterscheiden sich voneinander offensichtlich nur 
durch voneinander abweichende Definitionen des Gegenstandes, fiir den Schutz 
begehrt wird und die fiir die Merkmale dieses Gegenstandes venwendete 
Terminologie. Somit sind die Anspruche nicht knapp gefaBt. Femer mangelt es 
den Anspriichen insgesamt an Klarheit, da es aufgrund der Vielzahl unabhangiger 
Anspruche schwierig, wenn nicht unmoglich ist, den Gegenstand des 
Schutzbegehrens zu ermittein, und damit Dritten die Feststellung des 
Schutzumfangs in unzumutbarer Weise erschwert wird. 
Aus diesem Grund erfullen die Anspruche 1 , 8, 9, 10, 1 1 , 12 und 14 nicht die 
Erfordernisse des Artikels 6 PCT. 
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3.1 Urn diesen Einwand auszuraumen, scheint es zweckmafBig, einen geanderten 
Satz Anspriiche einzureichen, in dem der Gegenstand mit einem einzigen 
unabhangigen Anspruch definiert wird, wobei noch abhangige Anspruche 
hinzukommen konnen, die die fakultativen Merkmale abdecken (Regel 6.4 PCT). 
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VERTRAG UBgi^DIE INTERNATIONALE ZUSAMyfi^ 
^ AU|[p:M GEBIET DES PATENTWES^P 

PCT 

INTERNATIONALER RECHERCHENBERICHT 

(Artikel 18 sowie Regein 43 und 44 PCT) 



MARBEIT 



Aktenzeichen des Anmelders Oder Anwalts 
S04 P03 WO 


WEITERES MItteilung uber die Ubermittlung des intemationalen 

Recherclienberichts (Formblatt PCT/ISA/220) sowie, soweit 
VORGEHEN zutreffend, nachstehender Punkt 5 


Internationales Aktenzeichen 
PCT/DE 00/02930 


Internationales Anmeldedatum 
(Tag/Monat/Jahr) 

28/08/2000 


(Fruhestes) Prioritatsdatum (Tag/Monat/Jahr) 

28/08/1999 


Anmelder 

SQUID AG 



Dieser Internationale Recherchenbericht wurde von der Intemationalen Recherchenbeh6rde erstelit und wird dem Anmelder gem^B 
Artikel 18 ubermittett. Eine Kopie wird denfi Intemationalen Buro ubermittelt. 

Dieser internationale Recherchenbericht umfaBt insgesamt _5 Blatter. 

[X] Daruber hinaus llegt ihm Jewells eIne Kopie der in diesem Bericht genannten Unterlagen zum Stand der Technik bei. 



1 . Grundlage des Berichts 

a. Hlnsichdich der Sprache 1st die internationale Recherche auf der Grundlage der intemationalen Anmeldung in der Sprache 
durchgefOhrt worden, in der sie eingereicht wurde, sofern unter diesem Punkt nichts anderes angegeben ist. 

I [ Die internationale Recherche ist auf der Grundlage einer bei der Behorde eingereichten Ubersetzung der intemationalen 
Anmeldung (Regel 23.1 b)) durchgefuhrt worden. 

b. Hinsichtlich der in der intemationalen Anmeldung offenbarten Nucleotid- und/oder Aminosauresequenz ist die Internationale 
Recherche auf der Grundlage des Sequenzprotokolls durchgefuhrt worden, das 

I I in der intemationalen Anmeldung in Schrifiicher Form enthalten ist. 

zusammen mit der intemationalen Anmeldung In computerlesbarer Form eingereicht worden ist. 
bei der BehSrde nachtrSglich in schriftllcher Form eingereicht worden ist. 
bei der BehOrde nachtrSglich in computerlesbarer Form eingereicht worden ist. 

Die Erkiarung, daB das nachtrdgltch eingerelchte schriftliche Sequenzprotokoll nicht uber den Offenbarungsgehalt der 
intemationalen Anmeldung im Anmeldezeitpunkt hlnausgeht, wurde vorgetegt. 

Die ErkiSrung, dafB die in computerlesbarer Form erfaBten Infonmationen dem schriftlichen Sequenzprotokoll entsprechen, 
wurde vorgelegt. 

Bestimmte Anspruche haben sich als nicht recherchierbar erwiesen (siehe Feld I). 
Mangelnde Einheitlichkeit der Erfindung (siehe Feld II). 



2. 
3. 



□ 
□ 

□ 
□ 

□ 

□ 
□ 



4. Hinsichtlich der Bezeichnung der Erfindung 

fX] wird der vom Anmelder eingerelchte Wortlaut genehmlgt. 
I I wurde der Wortlaut von der BehOrde wie folgt festgesetzt: 



Hinsichtlich der Zusammenfassung 

I I wird der vom Anmelder eingerelchte Wortlaut genehmigt. 

wurde der Wortlaut nach Regel 38.2b) in der in Feld III angegebenen Fassung von der Beh6rde festgesetzt. Der 
I A I Anmelder kann der Behdrde Innerhalb eines Monats nach dem Datum der Absendung dieses intemationalen 

Recherchenberichts eine Stellungnahme vorlegen. 

Folgende Abbildung der Zeichnungen ist mit der Zusammenfassung zu verOffentlichen: Abb. Nr. _1 



I I wie vom Anmelder vorgeschlagen keine der Abb, 

Pn well der Anmelder selbst keine Abbildung vorgeschlagen hat. 
I I well diese Abbildung die Erfindung besser kennzeichnet. 
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INTERNATIONAl 



RECHERCHENBERICHT 



E 00/ 02930 



Feld III 



WORTLAUT DER ZUSAMMENFASSUNG (Fortsetzung von Punk! 5 auf Blatt 1) 



Computerbasieres Verfahren zur automatischen Aufbereitung von Daten 
biomagnetischer Felder, insbesondere von magnetokardiographischen Daten. 
Der Erfindung liegt die Aufgabe zugrunde, ein Verfahren zum Aufbereiten der 
Daten des erfa&ten biomagnetischen Feldes anzugeben, mittels welchem die das 
biomagnetische Feld hervorruf enden Strome durch Stromlinien in der Ebene derart 
visual isiert werden konnen, dafi sie einer visuellen Auswertung durch 
geschultes Personal oder einen Arzt in besonders einfacher Weise zuganglich 
werden . 
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EHENEl 



Intern ation ales Aktenzeichen 

_PCTj||^^0/02930_ 



Kategorie** 



Bezelchnung der Veroffentlichung, soweit erforderllch unter Angabe der in Betracht Komnnenden Teile 



TSUKADA K ET AL: "Noninvasive 
visualization of activated regions and 
current flow in the heart by analyzing 
vector components of a cardiac magnetic 
field" 

COMPUTERS IN CARDIOLOGY 1999. VOL.26 (CAT. 
N0.99CH37004), COMPUTERS IN CARDIOLOGY 
1999. VOL. 26, HANNOVER, GERMANY, 26-29 
SEPT. 1999, 

Seiten 403-406, XP002174332 
1999, Piscataway, NJ, USA, IEEE, USA 
ISBN: 0-7803-5614-4 

Seite 403, linke Spalte, Zeile 20 -rechte 
Spalte, Zeile 21 

Seite 404, linke Spalte, Zeile 4 -Seite 
406, linke Spalte, Zeile 19 

L.STADNYUK ET AL: "Interdipendence of 
Dynamics of the Echocardiography and 
Magnetocardiography Characteristics after 
Transesophageal El ectrophysl ol ogi cal 
Testing" 

REPORT ON IITH INTERNATIONAL CONFERENCE ON 
BIOMAGNETISM, 

1998, Seiten 136-140, XP002174333 
Sendai , JP 

Seite 136, rechte Spalte, Zeile 15 -Seite 
137, rechte Spalte, Zeile 36 

I.CHAIKOVSKY ET AL: "Ventricular 
Repolarization Disturbance Diagnostics in 
Chronic Ischemia Patients Derived from 
MCG" 

REPORT ON THE lOTH INTERNATIONAL 
CONFERENCE ON BIOMAGNETISM, 

1996, Seiten 210-213, XP002174334 
Santa Fe, USA 

Seite 210, linke Spalte, Zeile 27 -Seite 
212, linke Spalte, Zeile 19 
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1-14 
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INTERNATIONALER RECHERCHENBERICHT 



A. KLASSIFIZIERUNG DES ANMELDUNGSGEGENSTAKDES 

IPK 7 G06F17/00 G06F19/00 G01R33/035 



SEGENST 
F19/0( 



Interna tiona les Aktenzeichen 

PCT^k 00/02930 



A61B5/04 



Nach der Internationalen Patentklassifikation (tPK) oder nach der nationalen Klassifikation und der IPK 



B. RECHERCHIERTE GEBIETE 



Recherchierter Mindestprufstoff (Klassifikationssystem und Klassifikationssymbole ) 

IPK 7 G06F GOIR A61B 



Recherchierte aber nrcht zum Mindestprufstoff gehdrende Veroffentlichungen, sowelt diese unter die rscherchierten Gebiete fallen 



Wahrend der internationalen Recherche konsultierte elektronlsche Datenbank (Name der Datenbank und evti. verwendete Suchbegriffe) 

INSPEC, WPI Data, EPO-Internal , PAJ, MEDLINE 



C. ALS WESENTLICH ANGESEHENE UNTERLAGEN 



Kategorie* Bezek:tinung der Veroffentltehung. sowelt erfordedicti unter Angabe der In Betracht komnnenden Telle 



Betr. Anspruch Nr. 



NENONEN J ET AL: "Magnetocardiographic 

functional localization using current 

multlpole models" 

IEEE TRANSACTIONS ON BIOMEDICAL 

ENGINEERING, JULY 1991, USA, 

Bd. 38. Nr. 7, Selten 648-657, 

XP002174328 

ISSN: 0018-9294 

Seite 649, linke Spalte, Zeile 12 -Selte 

650, rechte Spalte, Zeile 9 

Seite 650, rechte Spalte, Zeile 29 -Seite 

651, linke Spalte, Zeile 41 

Seite 652, rechte Spalte, Zeile 7 - Zeile 
26 
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1-14 



Wettere Veroffentlichungen sind der Fortsetzung von Feld C zu 

entnehmen 
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Slehe Anhang Patentfamilie 



** Besondere Kategorien von angegebenen Veroffentlichungen 

■A' Veroffentllchung. die den allgemeinen Stand derTechnikdeflnlert, 
aber nicht als besonders bedeutsam anzusehen 1st 

'E' afteres Dokumenl. das jedoch erst am oder nach dem internationalen 
Anmeldedatum veroffentlicht worden ist 

•L* Veroffentlichung, die geeignet ist, einen Prioritatsanspruch zweifelhaft er- 
scheinen zu lessen, oder durch die das Veroffentlichungsdatum einer 
anderen im Recherchenbericht genannten Verdffentlk:hung belegt werden 
son Oder die aus einem anderen besonderen Grund angegeben ist (wie 
ausgefuhrl) 

'O' Veroffentlichung, die sich auf eine miindliche Offenbarung, 

eine Benutzung, eine Ausstellung oder andere Maf3nahmen bezieht 

•P' Veroffentlichung, die vor dem internationalen Anmeldedatunn, at)ernach 
denn beanspruchten Prioritatsdatum veroffentlicht worden ist 



"T' Spatere Veroffentlichung, die nach dem Internationalen Anmeldedatum 
Oder dem Prioritatsdatum veroffentlicht worden ist und mit der 
Anmeldung nteht kollidiert, sondem nur zum Verstandnis des der 
Erflndung zugrundeliegenden Prinzips oder der ihr zugrundeliegenden 
Theorie angegeben ist 

•X' Veroffentlichung von besonderer Bedeutung; die beanspruchte Erflndung 
kann allein aufgrund dieser Veroffentlichung nicht als neu Oder auf 
erfinderischerTatigkelt beruhend bet rachtet werden 

•Y' Veroffentlichung von besonderer Bedeutung; die beanspruchte Erflndung 
kann nicht als auf erfinde rise her TStigkeit beruhend betrachtet 
werden, wenn die Veroffentlichung mit einer oder mehreren anderen 
Veroffentlichungen dieser Kategorie in Verbindung gebracht wird und 
diese Verbindung fur einen Fachmann naheliegend ist 

Verdffentlichung. die Mitglied derselben Patentfamilie Ist 



Datum des Abschlusses der internatbnalen Recherche 
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Absendedatum des lntematk>nalen Recherchenberichts 



27/08/2001 
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Fax: (+31-70) 340-3016 
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14, rechte Spalte, Zeile 6 

RAMON C ET AL: "MCG simulations with a 
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